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Abstract 


This report describes work performed by United Technologies Research Center (UTRC) for 
NASA Langley Research Center (LaRC) under Contract NNL11AA06C. The objective of 
this program is to develop technology to reduce helicopter interior noise resulting from mul- 
tiple gear meshing frequencies. A novel active vibration control approach called Minimum 
Actuation Power (MAP) is developed. MAP is an optimal control strategy that minimizes 
the total input power into a structure by monitoring and varying the input power of con- 
trolling sources. MAP control was implemented without explicit knowledge of the phasing 
and magnitude of the excitation sources by driving the real part of the input power from 
the controlling sources to zero. It is shown that this occurs when the total mechanical input 
power from the excitation and controlling sources is a minimum. MAP theory is developed 
for multiple excitation sources with arbitrary relative phasing for single or multiple discrete 
frequencies and controlled by a single or multiple controlling sources. Simulations and ex- 
perimental results demonstrate the feasibility of MAP for structural vibration reduction of 
a realistic rotorcraft interior structure. MAP control resulted in significant average global 
vibration reduction of a single frequency and multiple frequency excitations with one con- 
trolling actuator. Simulations also demonstrate the potential effectiveness of the observed 
vibration reductions on interior radiated noise. 
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Chapter 1 

Executive Summary 


The objective of this work was to develop and demonstrate an approach for active vibration 
reduction on a rotorcraft structure to reduce noise in the passenger cabin of the rotorcraft. 
Vibration reduction at frequencies where gear-mesh tones and harmonics occur (500-2000 
Hz) due to structure-borne excitation from the main gearbox is of interest to the rotorcraft 
cabin noise community. The proposed work developed and demonstrated a control method- 
ology that is simpler to implement than traditional control approaches and potentially more 
robust and reliable for on-board implementation. Successful completion of this 3 year pro- 
gram will assist NASA in addressing the goals of the Quiet Cabin Technical Challenge within 
the Subsonic Rotary Wing project, where methods and concepts are sought-after that can 
be applied to or integrated into the structure to reduce rotorcraft cabin noise. 

Conventional active vibration control (AVC)/active noise control (ANC) typically incorpo- 
rates a sensor array for monitoring the vibratory or acoustic response, a controller, and 
control actuators. In many cases, applications of active control are hampered by the added 
system complexity, reduced reliability, and potential weight penalty. A novel approach to 
AVC/ANC, Minimum Actuation Power (MAP), is demonstrated that results in control of 
a single scalar objective, simpler and faster real time implementation through monitoring 
and adjusting the electrical input power of the control actuator and reduced hardware for 
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distributed vibratory response sensing, as compared with the conventional approach. MAP 
control features the following principle characteristics; 1) formulating the AVC process as 
real power control between vibration sources and actuators, 2) suppressing vibration by min- 
imizing the total input power as the sole objective function, and 3) eliminating the need for 
error sensing by using the reactive nature of actuator for self-sensing of vibratory response. 

The theory of Minimum Actuation Power (MAP) was first developed for a structure ex- 
cited by a single primary (disturbance) source and controlled by a single secondary (control) 
source. It was shown that the real input (mechanical) power in the secondary source goes 
exactly to zero when the total input power from both the primary and secondary source is 
minimized. This leads to a novel way to implement MAP control without explicit knowl- 
edge about the magnitude and location of the primary excitation. In order to validate the 
theoretical findings, a model for a simply supported plate excited by a primary piezo-electric 
(PZT) patch actuator and controlled by a secondary piezo-electric patch actuator was devel- 
oped. An electromechanical coupling model of piezoelectric patch actuator with a structure 
is developed and integrated with plate MAP model. A simple relation between the real 
electrical power and the real input power of the piezo-electric actuator was developed. This 
enabled implementing MAP control using the real electric power in the secondary actuator. 
The effectiveness of MAP control was demonstrated for resonant and off-resonant excitation 
frequencies. It was shown through simulations that optimizing the location of the secondary 
actuator can help control otherwise uncontrolled modes very effectively. The MAP theory 
was extended to multi-frequency excitation and simulations were presented to validate the 
theory. A test setup that mimics simply supported boundary conditions was built to vali- 
date the theoretical findings. To implement MAP using electrical power, tests were done to 
characterize the power consumption of the plate and the dielectric loss of the piezo-electric 
patch actuators. The MAP control algorithm was tested for on-resonance and off-resonance 
frequencies and effectiveness of MAP for vibration reduction. We demonstrated up to 47 dB 
vibration reduction at resonance and 7 dB vibration reduction at off-resonance, as measured 
by a single accelerometer at the plate center. There was good qualitative agreement between 
the experimental results and theoretical findings. 

The MAP control approach was extended to multiple primary excitations with zero relative 
phasing and multiple secondary excitations. It was shown that when the primary excitations 
have the same phasing, then the real secondary power is zero when the total input power is 
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minimized. However, when the primary excitations have arbitrary phasing, the secondary 
real power is not always zero when the total input power is minimized. One way the MAP 
condition is satisfied is when the secondary sources are located symmetrically with respect to 
the primary sources. In general, such a symmetry constraint is difficult to satisfy for arbitrary 
locations of the primary sources. An electro-mechanical model for a simply supported plate 
excited by multiple primary and multiple secondary PZT patch actuators was developed. 
Performance of MAP control as a function of the location and number of the secondary 
sources was assessed. It was seen that increasing the number of secondary sources improves 
vibration reduction. Also, multiple secondary sources can be useful when vibration reduction 
over a frequency band is sought. Experimental validation of MAP control was demonstrated 
in the lab with two primary and two secondary excitations using an aluminum plate as well 
as a composite panel. For example, we demonstrated 20 dB vibration reduction for the 
aluminum plate and 12 dB vibration reduction for the composite plate, as measured by a 
single accelerometer at the plate center. 

The feasibility of MAP control approach for reducing airframe vibrations due to the rotor- 
craft main rotor transmission and the reduction of corresponding sound levels in the passen- 
ger cabin was assessed using a more meaningful metric for vibration reduction, the global 
kinetic energy of the plate structure. MAP feasibility was assessed for multiple primary 
excitation sources with arbitrary relative phasing and with multiple frequencies controlled 
by multiple secondary sources. In order to circumvent the problem of placing the secondary 
actuators symmetrically with respect to the primary sources to satisfy the MAP condition, 
it was shown that placing the the secondary sources at the intersection of the node lines 
of the neighboring modes can guarantee that the secondary real power is zero at optimal 
conditions even when the primary sources have arbitrary phasing. An algorithm to compute 
the optimal secondary phase and magnitude was developed for resonant and off-resonant 
excitation frequencies. A stable iterative solution for finding the MAP solution for multiple 
secondary sources was developed. The stable iterative solution, however, is not feasible for a 
large number of secondary sources. An alternative feasible approach was developed to find 
the MAP solution when the secondary sources are large in number. A test rig representa- 
tive of a rotorcraft interior structure was designed, modeled and built. MAP control was 
successfully demonstrated with primary excitation using shaker, in addition to PZT, with 
multiple frequencies and multiple secondary actuators for control. It was shown that MAP 
control achieved significant vibration reduction with average global vibration reduction of 
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the order of 15-20 dB for single frequency with one secondary actuator. Also, a single sec- 
ondary actuator was optimized to control multiple frequencies, in some cases up to 3 with 
8-10 dB average global vibration reduction for each frequency. 

A fully coupled vibro-acoustic finite element model of a larger, more complex and modally 
denser box type structure was used as a quick evaluation of the potential of MAP to reduce 
noise in the passenger cabin due to the main rotor transmission excitation. It was shown 
that MAP control resulted in 6.1 dB reduction in the total kinetic energy of the structure 
while the averaged sound pressure levels (SPL) in the acoustic cavity dropped by 13.8 dB 
at resonant excitation. This preliminary analysis clearly shows the potential of MAP in 
reducing interior noise in a helicopter cabin. 
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Chapter 2 


Introduction 


The objective of the proposed work is to demonstrate active vibration reduction for interior 
noise reduction through the development of a control methodology that is simpler to imple- 
ment, more robust and reliable for on-board embedment, and requires less and/or low cost 
hardware. The focus is on demonstrating significant vibration reduction for mid-frequency 
(500 - 2000 Hz) cabin noise resulting from the gear-box transmission of structure-borne 
excitation into helicopter interiors. Active vibration/noise control for multiple tone struc- 
tural noise reduction conventionally employs a sensor array for monitoring the vibratory or 
acoustic response, a controller with a computationally expensive algorithm, and multiple ex- 
citation actuators. In many cases, applications of active control are hampered by the added 
system complexity, reduced reliability and weight penalty, especially for rotorcraft applica- 
tions where weight restrictions and reliability are stringent. The proposed active vibration 
and noise technology using minimum actuation power can significantly reduce the control 
hardware requirement and greatly simplifies the control algorithm by using a simple and 
unified objective control function and actuation self-sensing. 
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2.1 Motivation 


Interior noise is an increasingly important discriminator in the commercial helicopter market, 
with acceptable noise levels traditionally being achieved passively, albeit with substantial 
weight penalties. Increasing performance demands (i.e. , longer range, higher payloads) have 
driven the pursuit of lighter-weight and higher-performance solutions, often compromising 
the acoustic comfort inside the cabin. The main sources of high interior noise levels in a 
helicopter are the main rotor, tail rotor, and main gearbox. The main and tail rotors generate 
low frequency airborne noise, while the main gearbox generates mostly structure-borne noise 
from 500 - 2000 Hz, see Fig. 2.1. This mid- frequency range is the focus of the proposed 
research. 

The main gearbox dynamic excitations, generated by meshing gear pairs, are a significant 
source of vibration and cabin noise in a helicopter. The gearbox excitations are transmitted 
through to the fuselage via rigid connections that do not appreciably attenuate vibratory 
energy. These mid to high frequency vibrations typically include discrete gear-meshing fre- 
quencies that excite the frames and skins of the helicopter. Subsequently, the frame and skins 
of the fuselage transmit the vibrations produced by the main gearbox and radiate sound into 
the cabin interior. Figure 2.2 illustrates the key components of helicopter interior noise 
sources, propagation paths, and radiation to the interior. Over the past several years, stiffer 
and lighter weight composite structures have been replacing helicopter structures previously 
constructed from metal. Specifically, the aircraft skin which was traditionally aluminum 
has been replaced by composite sandwich panels which have a high Young’s Modulus and 
high flexural wave speeds that cause the panels to be efficient sound radiators in the fre- 
quency range of the main gearbox as shown in Fig. 2.3. For such composite panels, typical 
constrained layer damping passive noise treatment are not effective, necessitating either a 
different panel structure or an active control system to reduce the radiated noise. 
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Figure 2.1: Helicopter interior noise profile. 



Figure 2.2: Key helicopter interior noise sources, components, and structures. 

2.2 Background 


Controlling vibration propagation in the airframe is assumed to lead to improvements in 
passenger comfort and a better environment for on-board electronics. Active Control of 
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Figure 2.3: Measured radiation index for aircraft skin panels. 

Structural Response (ACSR) is a technology developed by Westland Helicopters Ltd., and 
places its actuation between two hard-points in the fuselage to produce equal and opposite 
forces between these two hard points. The ACSR approach utilizes these dual-point actua- 
tors to minimize the response of the fuselage to the non-propulsion vibratory loads. With a 
strut-mounted gearbox, and by adding magnetostrictive actuators, Sutton et al. reduced the 
transmitted kinetic energy to the fuselage of up to 40 dB at frequencies below 1250 [1]. Re- 
cently, this control of energy transmission through the strut was revisited with an active and 
periodic strut design [2] or refined control algorithms [3]. Unlike the Westland helicopters 
EH101, Sikorsky helicopters, such as the S-92A and UH-60, have the main gearbox bolted 
directly to the airframe instead of using struts. As a consequence, Sikorsky has developed 
an Active Vibration Control (AVC) architecture that involves the use of multiple airframe- 
mounted force actuators to provide canceling vibration, an approach that was successfully 
demonstrated in-flight on a Sikorsky UH-60L aircraft to cancel 4P vibration [4], Accelerom- 
eters mounted in the cockpit and cabin provided feedback to a controller to determine the 
actuator commands required to minimize vibration. A schematic of the configuration is pre- 
sented in Fig. 2.4. A similar AVC system is included as a baseline on all S-92A aircraft [5] 
and was selected for incorporation into the Army UH-60M, the Navy MH-60R and MH-60S 
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aircrafts [6, 7]. Finally, AVC is currently being incorporated into the S-76D helicopter and 
is also one of the enabling technologies for the Sikorsky X2 Technology Demonstrator. 



Feedback control architecture to produce vibratory loads which counteract those 
produced by the rotor system. 


Figure 2.4: Schematic of the AVC architecture developed by Sikorsky. 

Active Noise Control (ANC) is an adequate technology to reduce noise in the low- to mid- 
frequency range [8, 9, 10]. Using this approach in aircraft, many studies reported significant 
noise reduction, with most of the difference in performance coming from the differences in 
implementation and control algorithms [11, 12, 13, 14, 15]. Sikorsky Aircraft and the United 
Technologies Research Center have also been actively engaged in the development and testing 
of active noise control and active vibration control systems for a number of decades. In the 
early 1990s, Sikorsky embarked on a multi-year IR & D program to develop a flight-worthy 
ANC system to actively cancel gear- mesh noise inside the cabin interior. The classical ANC 
application with speakers as actuators proved inadequate for a typical helicopter cabin due 
to the high frequencies (> 700Hz) and corresponding large number of participating acoustic 
modes that would require that the number of speakers be at least equal to the number of 
relevant acoustic modes. The number of acoustic modes as a function of frequency for a 
cabin the size of S-76 helicopter is shown in Fig. 2.5 

Sikorsky instead pursued the more practical approach of using a pseudo choke-point method- 
ology to prevent the structure-borne gear-mesh vibratory energy from entering the cabin [16]. 
A schematic of this approach is presented in Fig. 2.6 and involved placing small resonant 
proof-mass actuators where the main transmission mounts to the airframe, and using micro- 
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Figure 2.5: Number of acoustic modes for a 5 ft x 6 ft x 9 ft enclosure. 


phones inside the cabin as error sensors. Two ANC flight tests were conducted on a Sikorsky 
S-76B commercial helicopter [17], the first of which is the first known successful flight test 
of a high frequency gear- mesh ANC system on a helicopter. The development flight tests 
demonstrated the capability of the gear- mesh ANC system to achieve tonal noise reductions 
of 18 dB on average and peak reduction by about 23 dB. Based on the success of these 
flight tests, a follow-on project was conducted to install a production- ready ANC system on 
a production S-76B helicopter with a VIP Cocoon Interior. The production demonstration 
flight test was successful in providing significant reductions in the primary gear-mesh tones. 

Finally, if the vibration propagation cannot be controlled efficiently or sufficiently, an alter- 
native strategy is to prevent the panels and skin from radiating noise inside the helicopter 
cabin. Active Structural Acoustic Control (ASAC) [18, 19] alters the noise radiation of the 
structure through the modification of the vibration pattern. Fuller et. al. [20] mounted 
piezoceramic actuators directly to an aircraft frame and with only two error microphones 
inside achieved a noise reduction of over lOdB inside the aircraft fuselage. However, it 
was shown that the radiated sound could be controlled without having error sensors in the 
acoustic field, having instead error sensors directly on the structure and a radiation model 
[21, 22]. As an implementation, the velocity feedback architecture was shown to introduce 
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Figure 2.6: Schematic of the Sikosrsky ANC architecture for helicopter interior noise reduc- 
tion. 

active damping by feeding the output of the sensor directly back to the piezoceramic patch 
actuator via a fixed gain [23]. Furthermore, when the sensor and actuator are collocated, the 
system is stable. This velocity feedback architecture allows for a decentralized configuration, 
and actuators of various shapes were successfully tested [23, 24], 

With passive noise treatments being too bulky or ineffective to control noise radiation from 
composite panels, active control of vibration and noise is a viable and attractive alternative. 
Until now, commercial deployment of active control in rotorcraft has been limited due to high 
costs, stability issues and overall system complexity. The proposed active control approach, 
to be described in sections to follow, addresses and avoids these limitations. 
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2.3 Program goals 


The objective of this work is to develop and demonstrate an approach for active vibration 
reduction on a rotorcraft structure that will ultimately reduce noise in the passenger cabin 
of the rotorcraft. The active control approach to be developed and demonstrated in this 
work is based on the concept of actuation power minimization. This is a decentralized 
control concept, meaning all sensing, control input calculation, and actuation is done local 
to a specific actuator attached to the structure. This local control approach reduces the 
need for extensive wiring to connect distributed actuators and sensors and eliminates the 
need for a centralized control processing unit that coordinates the efforts of many actuators. 
Although the simplicity of a decentralized control approach is desirable, the approach is not 
without challenges. These include ensuring that locally optimized control inputs produce a 
net reduction in a relevant global performance metric, such as sound power in the passenger 
cabin. A second challenge is to keep the actuator force requirements within realizable limits 
at the frequencies of interest. 

The three-year program is aimed at achieving the following specific overall objectives: 


1. Establish a fundamental understanding and expertise of the minimal power control 
methodology and demonstrate feasibility on a proof-of-concept test in Year 1. 

2. Develop modeling capability and a real time control algorithm and demonstrate the 
technology’s effectiveness and ease of implementation on a rotorcraft interior panel 
structure in Year 2. 

3. Demonstrate a technical path toward applying the methodology to the gear meshing 
noise of rotorcraft main rotor gear box and analytically prove the effectiveness in Year 
3. 


Specific year-wise goals as stated in the contract Statement of Work (SOW) are: 
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Phase 1 Theoretical Development and Proof -of -Concept Demonstration 


1.1 Develop the theory of minimal actuation power (MAP) vibration control of a single 
mechanical excitation (the primary excitation) of a structure using a single control 
actuator (the secondary actuator). 

1.2 Develop electromechanical models of the coupling between reactive secondary actu- 
ators and a host structure. 

1.3 Develop a real-time control algorithm to implement MAP for a single secondary 
actuator within the frequency range of 500-2000 Hz. 

1.4 Experimentally demonstrate MAP for a single primary excitation and single sec- 
ondary actuator on a simple structure at multiple frequencies. 

1.5 Present theoretical and experimental studies of single excitation/single control MAP 
to quantify control performance using metrics that include, but not be limited to, 
the reduction of a response variable, the secondary actuator force requirements, and 
the total input power from the primary and secondary excitations. 


Phase 2 Minimal Power Control for Multiple Actuators and Multiple Excita- 
tions 

2.1 Develop a model of MAP control for multiple primary excitations and multiple 
secondary actuators. This shall include consideration of system performance as a 
function of the phasing between primary excitations at a given frequency. 

2.2 Develop models of the electromechanical coupling of multiple reactive actua- 
tor/sensors to a host structure. 

2.3 Develop system simulation models of the MAP control approach with multiple reac- 
tive actuators to enable design and optimization of a multiple-actuator MAP control 
system. 

2.4 Experimentally demonstrate MAP vibration control on a rotorcraft interior panel 
structure with multiple primary excitations and multiple secondary actuators. 
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Phase 3 Feasibility Assessment of Applying MAP for Vibration Reduction of 
Airframe Structure 


3.1 Assess the feasibility of the MAP control approach for reducing airframe vibrations 
due to the rotorcraft main rotor transmission and the reduction of corresponding 
sound levels in the passenger cabin. This assessment shall include, but not be limited 
to, estimates of system performance versus secondary actuator force and location 
within the frequency range of 500-2000 Hz. 

3.2 Develop a power model of a generic reactive actuator for rigid body excitation. 

3.3 Experimentally demonstrate MAP control for vibration reduction on a mock-up of 
a rotorcraft airframe subject to simulated main rotor transmission excitation. 

3.4 Document the technical path to apply MAP control to reduce noise in the passenger 
cabin due to the main rotor transmission excitation. 


2.4 Proposed Active Vibration Control Approach 


Our active vibration control approach aims at reducing the interior noise through active 
vibration control of the helicopter’s fuselage and interior structure over the vibration fre- 
quencies that are efficient sound radiators. A unique and unconventional concept to perform 
the structural active vibration and noise control features the following principal characteris- 
tics: 


• Formulating AVC process as an optimal control strategy that minimizes the total 
input power into the structure from all the sources of excitation that includes both the 
primary and secondary sources. 

• Minimizing the total input power using the input power in the secondary sources only 
without any knowledge of the primary sources. 

• Eliminating the need for error sensor by using the reactive nature of actuator for self- 
sensing of vibratory response. 
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MAP active vibration control has the following advantages: 


• Simple implementation. A single scalar quantity, the real power in the secondary 
source, is used to minimize the total input power. 

• Global vibration reduction. The total input power is an indicator of the global vibration 
level. 

• Reduced hardware. Since sensing and actuation is done at the secondary sources, the 
approach does not need additional sensors. 

• Unconditionally stable. 


A comparison with a conventional configuration is shown in Fig. 2.7. It is envisioned that 
implementation into a helicopter interior would potentially consistent of active vibration 
sensing and control on skins and or frames resulting in lowered measured interior noise 
levels, see illustration in Fig. 2.8. Note that the architecture in Fig. 2.8 tries to reduce 
interior noise levels by directly attenuating the panel vibration whereas the architecture in 
Fig. 2.6 reduces interior noise by sensing the interior noise using an array of micro-phones 
and using that information to suppress the vibration at the source locations. 


2.5 Report Outline 


In §3, we develop the theory of MAP for a structure excited by a single primary excitation 
and controlled by a single secondary excitation. We develop an electro-mechanical model 
of a simply supported plate excited by a single primary PZT actuator and controlled by 
a single secondary PZT actuator, then extend MAP theory for multi-frequency excitation. 
Simulation and experimental results are presented to validate the theoretical findings at 
resonance and off-resonance conditions. 


39 



voliage/freq 


(a) 



(b) 

Figure 2.7: Comparison between (a) Conventional feedback control and (b) proposed MAP 
control. 

In §4, we extend the theory of MAP for a structure excited by a multiple primary excitation 
and controlled by multiple secondary excitations. We develop an electro-mechanical model 
of a simply supported plate excited by multiple primary and secondary PZT actuators. The 
effect of relative phasing between the primary sources on MAP condition is analyzed and 
the performance of MAP control as a function of the number and location of the secondary 
actuators is studied. Simulation and experimental results are presented to validate the 
theoretical findings at resonance and off-resonance conditions. 

In §5, we assess feasibility of MAP control for rotorcraft vibration reduction with realistic 
excitation. Specifically, we look at feasibility of MAP for multiple primary sources with 
arbitrary relative phasing controlled by multiple secondary sources. A test-rig that mimics 
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Figure 2.8: Schematic of the potential architecture for gear-mesh structural AVC for interior 
noise reduction. 

rotorcraft construction was designed, modeled and built. Experimental results are presented 
that validate MAP control with multi-frequency shaker excitation and controlled by multiple 
secondary actuators for resonance and off-resonance excitations. 

In §6, we present a FEM analysis on a simple box type interior structure and demonstrate 
the potential of MAP in reducing the interior noise by reducing the structural vibration. 

In §7, we summarize the key accomplishments and provide next steps. 
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Chapter 3 

MAP Control for Single Primary and 
Single Secondary 


In this chapter, we develop the theory of MAP for a structure excited by a single primary 
source and controlled by a single secondary source. In §3.1, we review the relevant literature 
in the areas of structural vibration control using power minimization. In §3.2, we discuss 
the general theory of input power minimization and present the key idea behind MAP. A 
model of simply supported plate excited by a primary PZT actuator and controlled by a 
secondary PZT actuator is developed. The MAP theory is developed for multi-frequency 
excitation and a special case of commensurate frequencies is considered. Simulations results 
are presented in §3.3 to validate the theoretical findings. The performance of MAP control is 
evaluated at resonant and off-resonant frequencies. In §3.4, a simple relation between the net 
electrical power in the secondary actuator and the net input (mechanical) power is derived 
that enables implementing MAP by monitoring the net electrical power. Experimental results 
are presented in §3.5 


42 


3.1 Literature Review 


This section summarizes the published works specific to our proposed goals of implementing 
MAP control strategy and using piezoelectric actuators for vibration control of a panel 
structure. 

The concept of power minimization for the sound radiation of multiple acoustic monopoles 
was presented by Nelson et al. [25] and it was found that in many instances, there exists 
a unique set of secondary sources of complex source strengths that minimizes the overall 
acoustic power in the system. Elliott et al. [26] further developed the methodology by 
introducing two control strategies, i.e. the minimization of total input power by primary 
and secondary actuators and the maximization of secondary actuator power absorption. For 
the case of single primary and single secondary acoustic sources, Johnson and Elliott [27] 
demonstrated that the real power of a secondary source was zero when optimally controlling 
a primary source in an enclosed field. In addition, these authors noted when the two sources 
were close to each other, greater global sound reduction could be achieved. As will be shown, 
the proposed MAP approach follows the same reasoning and expends on it and applies it to 
vibration control. 

Bardou et al. [28] applied a similar strategy of considering the total power that enters a 
structure to control the flexural vibration of plates. Bardou et al. showed that minimizing 
the total power in the plate yields greater overall power attenuation than what results from 
maximizing the secondary power absorption, for both an infinite and finite plate. Further- 
more, the greatest attenuation was achieved if the secondary sources were relatively close 
to the primary source, i.e. when the two excitation sources are close with respect to the 
wavelength of the flexural wave, and in the case of infinite plate, the primary and secondary 
sources cannot be too distant if global control is to be achieved. Benassi and Elliott [29] also 
noted the effect of distance on global control, and used idealized passive treatments with 
specified impedances as secondary reactive sources to attenuate the vibratory energy in a 
plate structure. 

Another major task in MAP development, the methodology of actuator/sensor interaction 
with the vibrating host structure and its power consumption has been broadly reviewed 
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and analyzed. Although an electromagnetic actuator, which is of reactive nature, is capable 
of performing the collocated sensing and actuation in principle, its low sensitivity, large 
weight and implementation issues disqualify it as the first choice for the MAP demonstration. 
The piezoelectric actuator, on the other hand, is considered the most promising device for 
demonstrating the principle of the MAP AVC because of its intimately close coupling to the 
structure and its distributed and non-invasive actuation. 

A lot of work on piezoelectric power consumption has been reported since the first publication 
of the piezoelectric collocated sensor actuator by Dosch [30]. Early research of piezoelectric 
power consumption as actuator consists of two distinct approaches described in Liang et 
al[31, 32], Zhou et al [33] and Hagood et al [34], These approaches examined the power 
consumption of actuators when used to excite the host structure and did not consider the 
case of active control. Liang developed an equation of the one-dimensional piezoelectric 
actuator. The response of the actuator was determined by prescribing boundary conditions 
to the actuator that coupled it with the host structure. The resulting explicit solution of 
the coupled electromechanical impedance has been widely used for structural health mon- 
itoring and actuation efficiency optimization in modal control and modal analysis. The 
other modeling approach developed by Hagood et al, established a relation of piezoelectric 
actuator mechanical outputs with the electric voltage and current inputs to the actuator. 
Both models are proven to be accurate. Another comprehensive study by Warkentin [35] 
presented a development of piezoelectric power consumption for active vibration control us- 
ing the analytical models developed by Hagood et al. The conclusions in these papers are 
that the electrical power of piezoelectric actuators is dependent on the mechanical motion 
of the structure and the electrical characteristics of the piezoelectric material and geometry. 
More recent work on the actuator power modeling includes [36, 37, 38]. Brennan et al. [36] 
developed an ANC model that shows that electric power of a piezoelectric actuator in active 
vibration control is independent of the mechanical motion of the host structure when the 
structure is under optimal control. While the work was intended to define the electric power 
requirement for the amplifier, the conclusion has significant physical implications relative 
to MAP for AVC. Namely, when complete control is reached, the structure is motionless, 
therefore the power consumption of the piezoelectric actuator is no longer a function of the 
mechanical motion but a function of the geometry, material properties of actuators, and the 
voltage and frequency of the control signal (power baseline). This principle sets the stage 
for the electric MAP baseline control which simply seeks a blocked boundary condition to 
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the actuator without the need for knowledge of the primary excitation. The current devel- 
opment will build upon that work and bridge the gap between the global vibration reduction 
and minimal electric power (baseline). The spatial optimization of the actuators for modal 
control is another important aspect to consider [39, 40]. 


3.2 Theory of Minimum Actuation Power (MAP) 


In this section, we develop the theory of minimum actuation power. In §3.2.1, we develop 
the theory of MAP by considering a general structure that is excited by a set of primary 
sources and controlled by a set of secondary sources. We consider point forces for sake 
of illustration. However, the sources can be either forces or moments or a combination of 
both. We derive the optimal secondary sources that minimize the total power input into the 
structure. For the case with one primary source and one secondary source, we show that 
at optimal condition the real power in the secondary source is exactly zero. In §3.2.2, we 
extend the theory of MAP for a simply supported plate excited by a primary PZT patch 
and controlled by a secondary PZT patch. 


3.2.1 Input power minimization 


Conside a two-dimensional structure (e.g., simply supported plate) that is excited by N 
normal point forces 

/ T (t) = [/i / 2 ••• f N ] (3-1) 

and let 


v T (t) = [vi V 2 ... V N ] 


(3.2) 
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be the the velocity vector consisting of the set of normal velocities at the point of application 
of the forces. The time averaged total input power due to all the forces is given as 1 

P = ^Re{F H V} (3.3) 

where F and V represent the complex amplitudes of f and v. The total input power can 
be expressed as 

P = e{F H MF} (3.4) 

where M is the mobility matrix. Using reciprocity, M is symmetric. In this case 

P = ^ F H Re{M}F = ^F h RF (3.5) 

where R is the real part of M. If the force vector consists of primary and secondary sources 


f t = {fTfJ) 


then R can be decomposed as 


R 


Rpp Rps 

R sp R'ss 


where R sp = Rl s Using (4.7) and (3.7), we can express the total input power as 


p = -\F R F + F R F + F R F 4- F R F 1 

1 r) L X p - Lg 'pp J - p I x p - Ll/ ps- M - S ' x s - L * J Sp- L p l x s J -*'SS- a - Sj 


(3.6) 

(3.7) 


(3.8) 


Note that P is a quadratic function of F s . Minimizing P with respect to F s , we get [41] 

(3.9) 


popt — D ID p 

-*■ S J -*'SS J -*'Sp- L p 


1 The averaged product of two sinusoidally varying quantities represented in complex form reduces to the 
real part of the product of their complex amplitudes 

1, 


ui (t)u 2 (t) = -R eft/rRJ} 
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Key idea for MAP control 


Let us consider the case where there is one primary and one secondary source. In this case, 
the optimal secondary source is given by 


popt _ 


Rsp T? 

JD ^ P 


The real power in the secondary source is given by 

Ps — -Re{F:M sp F p + F*M SS F S } 

Using (3.10), the real power in the secondary source at optimal value is 


R 


Rl 


P: pt = -Re{-^F;M sp F p + ^F;M ss F p } 


1, R, 


R 2 


2 v p 1 p 1 L sp ± p w ± p 1 L ss ± pi 


= 0 


We make the following important observation: 


(3.10) 


(3.11) 


(3.12) 


For one primary force controlled by one secondary force, the real power P s in the secondary 
source is zero when F s = Ff pt where Ff pt is the optimal value that minimizes the total input 
power P 


This suggests a way to control the structure by only tuning the power in the secondary 
actuator without explicit knowledge about the primary source. The control strategy is to 
drive the real power in the secondary actuator to zero. This is the underlying principle of 
MAP control. 
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3.2.2 MAP theory for a simply supported plate actuated by two 
piezoelectric patches 

We investigate active control of the transverse vibrations of a simply supported plate using 
the Minimum Actuation Power (MAP) technique. The plate is excited by a primary PZT 
patch actuator and controlled by a secondary PZT patch actuator. Analytical expressions 
for the total input power due to the primary and secondary patches in terms of the equivalent 
patch angular mobilities are presented. The optimal value of the secondary moment that 
minimizes the total input power is computed. We show that at optimal value of the secondary 
moment, the total input power due the secondary actuator is zero. This key observation is 
the basis of MAP control. 


System description 


Consider a simply supported rectangular plate that is actuated by two identical PZT patch 
actuators as shown in Fig 4.1(a). Each actuator is driven by a voltage source. The primary 
patch acts as the disturbance while the secondary patch is the control actuator. Let L x , L y 
and H be the length, width and thickness of the plate respectively. The plate is assumed 
to be isotropic. The PZT patches are bonded to the plate surface and are assumed to be 
square with width a and thickness h. When an electric field is applied to the patch in the 
z direction, the patch expands and contracts in both the x and y directions. The resulting 
strain field generates a distributed moment along the edges of the patch, labeled 1 to 4, as 
shown in Fig 4.1(b). We assume that the patch is isotropic. Also, the patch is assumed to 
be thin and the patch inertia and stiffness are ignored in the analysis. 


Total input power from primary and secondary actuators 


Let M p and M s be the moments per unit length (complex amplitude) generated by the 
primary and secondary patch respectively. Let the primary and secondary patch edges be 
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Figure 3.1: (a) Schematic of a simply supported rectangular plate excited by two PZT 
actuators (b) Distributed moment generated by a patch. 

discretized into an equal number of elements and let Ae be the length of the edge element. 
Then, the primary and secondary moments generated by each element are 

M p = Ae M p 

M s = A eM s (3.13) 

The total power input from the primary and secondary patches can be written as 

P = P p + P s (3-14) 

where P p and P s are the total power inputs from the primary and secondary patches respec- 
tively. 

Let us assume that the primary and secondary patch edges are discretized into an equal 
number of elements, say q. Let Pi and s* represent the coordinates of the center of the i th 
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primary and secondary patch element respectively where i — 1, . . . , q. Let M p and M s be the 
primary and secondary moments acting on each element respectively. Consider the point Pi 
on the primary patch. The angular velocity at this point is a sum of the angular velocity 
induced by the primary and secondary patch moments and can be expressed as 

w Pi = e w kps + E w -“. < 3 ' 15 ) 

3 3 

where W PiPj is the angular velocity induced at Pi due to the moment M p at location pj and 
W PiS . is the angular velocity induced at pi due to the moment M s at location Sj. The angular 
velocities W PiPj and W PiPj can be expressed as 


W PiPj — PpiPj 

WpiSj — P pi SjM s 


where P PiPj and P PiSj are the cross mobility terms given by [28] 


Pp iPj = (wul + 


4 ( Pi ) 3pn j ( Pj 


P, 


PiSj 


m [(c^ - c o 2 ) 2 + uPr} 2 \ 

ST ( 2,-/2 2^ 4 {Pi) Ipn* (sj) 

> (r)UJUJ n + JU (Uj n - u )) — — — — 

^ rn < - UJ 2 2 + u^p 2 } 


72—1 

OO 


(3.16) 


(3.17) 


where p is the structural damping, m is the plate mass, c o is the excitation frequency (assumed 
same for primary and secondary excitation) and cu n is the 77 th resonant frequency given by 


^ n 


EH 2 


12(1 - u?)p 




(3.18) 


where n refers to the mode index ( 77 , 1 , 77 . 2 ), E is the modulus of elasticity, /1 is the Poison 
ratio and p be the density of the plate material. The function ijjn {p) is given by 
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where (x, y ) are the coordinates of the point p and 5 p is the angle made by the moment axis 
with the x axis. Note that for the primary (or secondary) patch, S Pi (or5 Si ) = 0, 7 t/2 , n and 
37 t/2 for all elements i along edges 1, 2, 3 and 4 respectively. The total power input due to 
the primary patch is given by 


Pr, 




(3.20) 


Using (A.l) and (A. 2) in (A. 6), the total input power due to the primary patch can be 
expressed as 

P p = ^Re ( m;p vp m p + m;p ps m s ) (3.21) 

where P pp represents the equivalent angular mobility due to the primary patch given by 


Ppp ~ PpiPj (3.22) 

i 3 

and P ps represent the cross angular mobility given by 

p,TE p w < 3 - 23 ) 

* 3 

Similarly, the total power input due to the secondary patch can be shown to be 

P s = ^ Re ( M* s P sp M p + M*P SS M S ) (3.24) 

where P sp = P ps represents the cross angular mobility and 

p»=EE p «i I 3 - 25 ) 

i 3 

represents the equivalent angular mobility due to the secondary patch. It can be shown that 


Pp — 2^ e (AfpPppMp + M*P ps M s ) 

(3.26) 

P, = ^Re (M*P sp M p + M* S P SS M S ) 

(3.27) 


51 


where P pp and P ss are the equivalent angular mobility terms due to the primary and secondary 
patch respectively while P ps and P sp are the cross angular mobility terms. By reciprocity, we 
have P ps = P sp . In other words, the total input power of a simply supported plate excited by 
two PZT patches is equivalent to the total input power of a simply supported plate excited 
by two point moments M p and M s with equivalent point and transfer mobilities P pp , P ss , P ps 
and P sp . The reader is referred to §B for a detailed derivation of the relations in (3.26) 
and (3.27). From (4.5), (3.26) and (3.27), the total input power from both the primary and 
secondary patch is given by 

P = e(M;P pp M p + M;P ps M s + M* s P sp M p + M* s P ss M s ) (3.28) 

Using reciprocity of the cross mobility terms, (4.6) can be further simplified as 

p = ^ (m;s pp m p + m;s ps m s + m*s sp m p + m*s ss m s ) (3.29) 

where S pp , S ps , S sp and S ss are the real parts of P pp , P psi P sp and P ss respectively. Note that 
S ps = S sp . 


Total input power minimization (MAP) control 


We seek to find the optimal secondary excitation that will minimize the net input power P. 
Equation (3.29) is a quadratic equation of the secondary moment M s and therefore a unique 
minimum exists. Minimizing the total input power P with respect to M s , we get the optimal 
secondary moment as a function of the primary moment and the angular mobility terms 

M°f = -^M p (3.30) 

^SS 

Since S sp and S ss are real, (3.30) shows that M° pt is either in phase or out of phase with M p 
depending on the sign of — S sp / S ss . An important implication of the above observation is that 
one only needs to vary the amplitude of the secondary moment to search for the optimal value. 
Equations (3.30) show that in order to compute the optimal secondary moment, we need 
information about the primary excitation. From an experimental point of view, frequency 
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of the primary excitation can be easily extracted. However, extracting the magnitude of the 
primary source is difficult. In that case, we cannot calculate the optimal values M° pt . In 
what follows, we show how the net power in the secondary actuator can be used as a control 
variable to implement MAP control without explicit knowledge about the strength of the 
primary excitation. Using (3.30) in (3.27), the real power in the secondary source at optimal 
value is given by 

P7 1 = + 

= \ + H 

= 0 (3.31) 

Thus, for a plate excited by a primary PZT actuator and controlled by a secondary PZT 
actuator, the real power P s in the secondary actuator is zero when M s = M° pt where M s opl is 
the optimal value of the secondary actuator moment that minimizes the total input power P 
from both the primary and secondary actuators. This suggests a way to control the structure 
by only tuning the power in the secondary actuator without explicit knowledge about the 
magnitude of the primary source. The control strategy is to drive the real power in the 
secondary actuator to zero. 


3.2.3 MAP control for multifrequency excitation 

We show that the theory of MAP derived for single frequency excitation also holds when 
you have multi-frequency excitation provided the frequencies are commensurate, i.e., integer 
multiples of a fundamental frequency. In section 3.2.3, we develop the theory of MAP 
for multi-frequency excitation by considering point forces acting on the structure. Though 
we consider point forces for sake of illustration, the results hold for distributed moment 
excitation as generated by a PZT patch. In section 3.3.7, we present simulation results for 
a simply supported plate excited by a primary PZT actuator and controlled by a secondary 
PZT actuator where the primary and secondary excitations are multi-frequency. 
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Let 

n 

/(*) = £/»(*) (3-32) 

k = 1 

be a point force excitation acting on the structure where f k is the k th component with 
frequency uj k - Let us assume that cu^’s are commensurate. Using superposition, the response 
at the excitation location can be expressed as 

n 

v (t) = ^2 v k{t) (3.33) 

k = 1 


where Vk is the response corresponding to the component fk of the excitation. The time 
averaged input power due to the excitation is given by 


P = f(t)v(t) = - 


[ (l>w) (x>w) dt 


(3.34) 


where T is the period. 


Since the frequencies are commensurate, Parseval’s lemma [42] gives 


n ^ pT n 

P = fk(t)v k (t ) dt = V P k (3.35) 

k = i 1 k = i 

where 

1 f T 

= — J fk(t)v k (t ) dt (3.36) 

is the time-averaged input power corresponding to oj k ■ Equation (3.35) will hold even when 
the system is excited by a primary and secondary force given by 2 


f P (t) = ^f P k(t) 

k = 1 
n 

f s {t) = ^ 3 - 37 ) 

k=\ 

2 In general, (3.35) holds if the structure is excited by a set of primary and secondary forces with similar 

harmonic components 
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where f pk and f sk are harmonic with frequency c o k . Note that in this case, P k denotes the 
total input power due to f pk and f sk . It follows from (3.35) that 

min P = {min P\ , min P 2l . . . , min P n } (3.38) 

That is, the total input power can be minimized by minimizing the input powers correspond- 
ing to each frequency. The optimal value of the secondary force that minimizes the total 
input power P is given by 

n 

/rw=E 4 tw- (3.39) 

k = 1 

where the optimal complex amplitude F° k l of (t) is given by 

FT = -R7,'MR, r MF£ V k (3.40) 

where F" k ' is the complex amplitude of f p k(t). Equation (3.40) follows from the MAP results 
for single frequency excitation. Thus, the MAP theory can be extended to multi-frequency 
excitation. 


Next, we show that minimization of total input power is equivalent to minimization of the 
net vibrational energy in the case of multi-frequency excitation. For example, consider a 
rectangular plate that is excited by primary and secondary forces as given by (3.37). By 
superposition, the velocity at any point (x, y ) on the plate can be expressed as 

n 

v(x, y,t)=^2 v k (x, y, t) (3.41) 

k= 1 


where v k (x, y : t ) has frequency oj k and is a contribution of both f pk (t ) 
vibrational energy of the plate is given by 


E ki n(t) = -m 


v(x, y, t ) 2 dx dy 


'0 Jo 

where m is the mass per unit area. The average vibrational energy is 


Ekin 


r L x rF r T ( n \ ( n 

Ilf \ 

Jo Jo Jo \ k=1 J \ k=l 


Vk(x,y,t ) 


and fskip)- The global 


(3.42) 

given by 

dtdxdy (3.43) 


55 


(3.44) 


Using Parseval’s lemma gives 



k = 1 




v k (x, V: t ) 2 dt dxdy = Y^ E kin 

k = 1 


where 


E k in = 


2 T 


y, t ) 2 dt dx dy 


(3.45) 


'0 ^0 ^0 

is the average vibrational energy corresponding to a;*,. Using the theory of MAP for single 
frequency, we have 

— 

min P k = min E kin (3.46) 

Using (3.38), (3.44) and (3.46), we get 

min P = min E^ m (3.47) 


Equation (3.47) shows that minimizing the average input power is equivalent to minimizing 
the average vibrational energy even when the excitation is multi-frequency provided that the 
frequencies are commensurate. 


3.3 Simulations 


We present simulation results to validate the analytical findings. The effectiveness of the 
control depends not only on the position of the secondary patch with respect to the mode 
shapes but also whether the excitation frequency is close to a resonant frequency or away 
from resonance. We present simulation results that help quantify the control effectiveness 
as a function of the secondary position and the excitation frequency. 


3.3.1 Simulation parameters 

A schematic of plate system for simulation studies is shown in Fig. 4.2. Note that ( x p ,y p ) 
and ( x s ,y s ) represent the coordinates of the centers of the primary and secondary patch 
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Figure 3.2: Schematic of the plate system for simulation. 


respectively. The location of the primary patch is fixed. We assume that the plate is made 
of aluminum. Table 4.1 lists geometrical and material properties for the plate and PZT 
patches. 


3.3.2 Vibration response under primary excitation 

We first study the response of the plate under primary excitation only. The secondary 
actuator does not apply any control. The coordinates of the center of the primary actuator 
are (x p ,y p ) = (0.275,0.225). The frequency range for the excitation frequency is chosen to 
be 100-2000 Hz. We assume that the voltage applied to the primary actuator is such that 
it generates a unit moment at each element, that is, M p = 1 Nm. There are 13 resonant 
frequencies in this frequency range: cu n (Hz) = 620, 724, 826, 976, 1039, 1222, 1395, 1413, 
1474, 1626, 1862, 1893 and 1982. Figure 3.3 shows the average vibrational energy of the 
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Table 3.1: Simulation Parameters 


Variable 

Value 

Lx 

0.38 m 

Ly 

0.35 m 

H 

4.2e-3 m 

E 

70e9 Pa 

V 

0.01 

P 

2700 kg/m 3 

P 

0.33 

a 

0.05 m 

h 

0.13e-3 m 

Y u 

6.6el0 Pa 

d'M 

-190e-12 m/V 

^33 

1.53e-8 F/m 

Up 

0.35 

5 

0.02 


plate when the primary excitation is applied. Interestingly, it can be seen from Fig. 3.3 
that the resonances at 1395 Hz and 1893 Hz are not excited. This is due to the fact that 
the location of the primary patch is close to the nodal lines of the mode shapes at 1395 Hz 
and 1893 Hz. Figure 3.4 shows the modes shapes at 1395 Hz and 1893 Hz. The position of 
the primary patch is also shown and it can be clearly seen that in both cases, the primary 
patch is located at the intersection of two nodal lines. Therefore, the primary excitation is 
not effective in exciting the modes at 1395 Hz and 1893 Hz. 


3.3.3 Vibration response with MAP control 

The center of the secondary actuator is located at ( x s ,y s ) = (0.075,0.075). Note that this 
location is not optimized for any frequency range. Figure 3.5 shows the average vibrational 
energy of the plate with and without MAP control. As seen, the control is most effective 
close to each resonance (the greatest difference in vibrational energy between control off and 
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control on). Note that the inodes at uj n = 1413 Hz and u j n = 1862 Hz are not controlled. 
This is because the location of the secondary actuator is not optimal for controlling these 
resonant modes. This is evident from Fig. 3.6 which shows the mode shapes for 1413 Hz 
and 1862 Hz, and at both resonances the secondary actuator is on the node lines and away 
from the mode shape peaks. 


3.3.4 MAP control at resonant and off-resonant frequencies 

The simulations show that the real part of the secondary power is zero when the total input 
power is minimum at both resonant and off-resonant frequencies. The key difference is that 
at resonant frequencies, the vibrational energy is also minimum when the total input power is 
minimum. However, at off-resonant frequencies, the vibrational energy is not truly minimum 
when the total input power is minimum (see Eq. B.ll). A more detailed analysis of the 
global vibrational energy and its relation to the total power input is presented in § B. 



Frequency (Hz) 

Figure 3.3: Vibrational energy of the plate under unit primary excitation. 
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to =1395 Hz Mode = 3,3 to =1893 Hz Mode = 4,3 

n n 



Figure 3.4: Resonance modes for cu n =1395 Hz and (n n =1893 Hz. The primary actuator 
location is shown by the black box. 



Figure 3.5: Vibrational energy of the plate with MAP control (dashed line) and without 
MAP control (solid line). 
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to =1413 Hz Mode = 1,4 
n 




Figure 3.6: Resonance modes for o; n =1413 Hz and u; n =1861 Hz. The secondary actuator 
location is shown by the black box. The primary actuator location is shown by the dashed 
box. 

For the sake of illustration, we focus on two frequency ranges: 550-750 Hz and 1500-1800 
Hz. Fig 3.7 shows a zoomed-in view of the vibration response with and without control 
for the frequency range 550-750 Hz. Two resonant frequencies (620 Hz and 724 Hz) and 
one off-resonant frequency (680 Hz) were chosen to study the relation between vibrational 
energy, total input power and real power in secondary actuator. At each frequency, we vary 
the secondary moment M s around the optimal value M° pt by defining 

M s = aM° pt (3.48) 

where a is a real scalar. Note that a = 0 implies zero control and a = 1 implies optimal 
control. 

Figures 3.8 and 3.9 show the variation of vibrational energy, total input power and real power 
in secondary actuator for the resonant frequencies 620 Hz and 724 Hz respectively. As a 
is varied from zero (no control), the vibrational energy and the total input power (primary 
plus secondary) decrease till a = 1 and then increase. The real secondary power decreases 
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650 

Frequency (Hz) 


Figure 3.7: Vibrational energy of the plate with MAP control (dashed line) and without 
MAP control (solid line) for the frequency range 550-750 Hz. Two resonant frequencies (620 
Hz and 724 Hz) and one off-resonant frequency (680 Hz) are marked by the blue triangles. 

and crosses zero at a = 1. In other words, at a = 1, both the vibrational energy and the 
net input power are minimum while the real power in the secondary is zero. On the other 
hand, Fig. 3.10 shows that when the excitation frequency is away from resonance, the total 
input power is minimum and the real power in the secondary actuator is zero when a = 1 
but the vibrational energy is not necessarily minimum at a = 1 though it is close to the true 
minimum. 

Similar conclusions can be drawn when we consider the frequency range of 1500-1800 Hz. 
Figure 3.11 show the vibration response with and without control for the above frequency 
range. We choose one resonant frequency (1626 Hz) and two off-resonant frequencies (1525 
Hz and 1750 Hz) for examining the variation of vibrational energy, total input power and 
the real secondary power as a is varied. 

Figures 3.12-3.14 show that at resonant frequencies, the real power in the secondary is zero 
when the total input power is minimum. Also, the vibrational energy is a minimum at a = 1. 
However, at off-resonant frequencies, the real power in the secondary is zero when the total 
input power is minimum but the vibrational energy is not necessarily minimum at a = 1. 
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Figure 3.8: Variation of the vibrational energy (top), the total input power (middle) and the 
real power in secondary actuator (bottom) with a for oj = 620 Hz. 


At off-resonance frequencies, setting the real power in the secondary actuator to zero will 
not necessarily produce the greatest attenuation in vibrational energy. 
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Figure 3.9: Variation of the vibrational energy (top), the total input power (middle) and the 
real power in secondary actuator (bottom) with a for oj = 724 Hz. 

3.3.5 Optimal secondary location for a frequency band 

Figure 3.5 shows that the secondary actuator position (x s , y s ) = (0.075, 0.075) is not optimal 
for the resonant modes at u> n = 1413 Hz and u n = 1862 Hz. If the secondary actuator location 
is chosen carefully, we can expect to control the above modes. We consider the frequency 
band 1300-1900 Hz to demonstrate optimization of the secondary actuator position. First, 
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Figure 3.10: Variation of the vibrational energy (top), the total input power (middle) and 
the real power in secondary actuator (bottom) with a for uj = 680 Hz. 


we note that the peaks and valleys of the mode shape for a given resonance mode are optimal 
locations for controlling that mode. Let <t>i{x,y ) denoted the mode shape for the resonant 
frequency cu* in the frequency band 1300-1900 Hz. We define the following optimization 
function that is a weighted sum of the mode shapes: 

$ = - ^ A abs(^(a;, y )) (3.49) 
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Figure 3.11: Vibrational energy of the plate with MAP control (dashed line) and without 
MAP control (solid line) for the frequency range 1500-1800 Hz. One resonant frequency 
(1626 Hz) and two off-resonant frequencies (1525 Hz and 1750 Hz) are marked by the blue 
triangles. 


where fa is some weight. The optimal locations for controlling the resonant modes in the 
frequency band are given by the minima of the function <f>. For the frequency band of interest, 
that is 1300-1900 Hz, there are resonances at = 1413 Hz, 1474 Hz, 1626 Hz and 1862 Hz. 
Figure 3.15 shows the plot for <3> where the minima are given by the dark blue areas. For this 
example, we weight the mode shapes equally, that is, fa — 1 for i — 1, . . . , 4. Two locations 
of the secondary actuator are shown. One of the locations, shown by the dashed box, is the 
non-optimal location given by ( x s ,y s ) = (0.075,0.075). The other location, shown by the 
solid box, is chosen to be one of the four minima and is given by ( x Sl y s ) = (0.125,0.135). 
Note that this optimal location is farthest from the primary location ( x p ,y p ) = (0.25,0.20). 
Figure 3.16 compares the vibration response when the secondary actuator is located at 
(x s ,y s ) = (0.075,0.075) and ( x s ,y s ) = (0.125,0.135). Figure 3.16 clearly shows that by 
optimizing the location of the secondary actuator, we are able to control all the modes in 
the given frequency band. 
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Figure 3.12: Variation of the vibrational energy (top), the total input power (middle) and 
the real power in secondary actuator (bottom) with a for uj = 1626 Hz. 

3.3.6 MAP control as a function of separation between primary 
and secondary patch 

The effectiveness of MAP control not only depends upon the location of the secondary 
actuator with respect to mode shapes, as we saw in the last section, but also on the distance 
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Figure 3.13: Variation of the vibrational energy (top), the total input power (middle) and 
the real power in secondary actuator (bottom) with a for uj = 1525 Hz. 


from the primary actuator. Let the reduction in vibrational energy be defined as: 


77 

~T^ re< ^ ^ kin ^ kin 

^ kin “ t=t 

^ kin 


(3.50) 


rpfj ~r==opt 

where E kin is the normalized reduction in the vibrational energy. E kin and Eki n are the 
vibrational energies of the plate with and without MAP control respectively. Note that 
E kin is a measure of the effectiveness of MAP control. Let the excitation frequency be 800 
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Figure 3.14: Variation of the vibrational energy (top), the total input power (middle) and 
the real power in secondary actuator (bottom) with a for uj = 1750 Hz. 


Hz (off- resonant). Figure 3.17 shows the general trend of E kin with distance between the 
primary and secondary actuator. In general, the effectiveness of MAP control decreases as 
the distance between the primary and secondary actuator increases. The control is most 
effective when the actuators are collocated. However, this is not feasible in practice and the 
best one can get is sub-optimal performance. Note that at nodal lines, the effectiveness of 
any control is zero as seen by the dips in Fig. 3.17. 
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Figure 3.15: Contour plot of the optimization function <3>. Dark blue areas indicate where 
the function is a minimum. The dashed black box shows the original secondary patch 
location (. x s ,y s ) = (0.075,0.075) while the solid black box shows the optimized location 
(. x s ,y s ) = (0.125,0.135). The primary actuator is shown by the solid brown box (top right). 

3.3.7 Multifrequency excitation 

We present results for MAP control when the primary excitation is multi-frequency. The 
primary patch location is given by (x p ,y p ) = (0.275,0.225) while the secondary patch is 
located at (x s ,y s ) = (0.075,0.075). Without loss of generality, assume that the primary 
excitation has two frequencies component uq = 1000 Hz and uj 2 — 1500 Hz. Note that uq 
and ou 2 are commensurate (integer multiples of some base frequency). Let M pl = M p2 = 1 
where M p \ and M p2 are the complex amplitudes of the primary moment corresponding to uq 
and uj 2 respectively. Let the components of the secondary moments be defined as a fraction 
of the corresponding optimal moments as: 

M s i = aM°f 

M s2 = aMf (3.51) 

where M^f' and are given by (3.30) with u replaced by ui and u 2 respectively. Recall 
that a is a real scalar where a = 0 implies no control while a = 1 implies optimal control. 
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Figure 3.16: Vibrational energy of the plate. No Control (solid); Secondary actuator at 
(. Xsifjs ) = (0.075,0.075) (red dashed line); Secondary actuator at (x s ,y s ) = (0.125,0.135) 
(blue dashed line). 



Figure 3.17: Vibrational energy reduction as function of the distance between the primary 
and secondary patch for c o = 800 Hz. 

The top plot in Figure 3.18 shows the variation of the total vibrational energy E with a. 
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Figure 3.18: Variation of the vibrational energy (top), the total input power (middle) and the 
real power in secondary actuator (bottom) for multi-frequency excitation with aq = 1000 
Hz and ca 2 = 1500 Hz. The bottom plot shows the components of the real power in the 
secondary actuator for oj\ = 1000 Hz (solid) and ca 2 = 1500 (dashed). 

Note that E = E\ + E 2 . The middle plot shows the variation of the total input power 
P = Pi + P 2 while the bottom plot shows the variation of P sl and Ps2 , components of the 
real power in the secondary actuator corresponding to uq and ca 2 respectively. The plots in 
Fig. 3.18 clearly show that when the total input power is minimum, the real power in the 
secondary actuator corresponding to each frequency is zero. Also, it can be seen that the 
total vibrational energy is minimum when the total input power is minimum. 
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3.4 MAP Implementation with Electric Actuator 


In an actual experiment, the net input power from the secondary actuator is hard to measure. 
Rather, one can easily measure the electrical power by monitoring the voltage and current 
in the patch. We develop the model for electromechanical coupling between the secondary 
actuator and the host structure. We relate the real electrical power in the secondary actuator 
to the net input power through the dielectric loss. We show how the real electrical power in 
the secondary actuator can be used to implement a real-time MAP control algorithm. 


3.4.1 Electromechanical coupling between the PZT actuator and 
the structure 

The two dimensional strain field developed in the PZT patch actuator is commonly ex- 
pressed by the following constitutive equation governing the electromechanical interaction 
of piezoelectric material: 
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(3.52) 


where e x and e x are the strains in the x and y directions respectively, a x and cr x are the 
stresses in the x and y directions respectively, E 3 is the electric field applied to the patch in 
the z direction, d 3 \ and d 3 2 are piezoelectric coupling coefficients, Yu is the elastic modulus 
and u p is the Poisson ratio of the piezoelectric material. Equation (3.53) can be expressed 
in terms of the stress field as 
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The electric displacement in the z direction is expressed as [43] 


(3.53) 


D 3 — s 33 E 3 + d 3 i<j x + d 3 2<J y 


(3.54) 
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where £33 is the complex dielectric constant of the piezoelectric material. Assuming the PZT 
is isotropic, that is d 3 1 = d 32 , we get 


D 3 =[e 33 -^-)E 3+ d J^( £x + s> ) 


J- I/p J A. Vp 

If V is the voltage applied to the PZT patch, we have 


(3.55) 


where h is the thickness of the PZT patch. Using (3.56) in (3.55), we get 

2 4iYn\ V , d 31 Y u 


D 3 — £33 


1 — v n 


- + 


1 — U r 


(E x + £y) 


The electric current in the PZT is given by 


n %2 ry 2 


I = juj 


D 3 dx dy 


J xi Jy 1 

where x 2 , y 1 and y 2 are the patch coordinates. Using (3.57) in (3.58), we get 
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where a is the width of the patch. The strain rates can be written as 

H + h d 2 v 

JW£ X = 


JUJ£y = 


2 dx 2 
H + h d 2 v 


2 dy 2 

where v is the transverse velocity of the plate. Using (3.60) in (3.59), we get 
2d 2 31 Y n \ a 2 V , d 31 Y n H + h f X2 f V2 f d 2 v , <9V 
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Consider the first term in the integral in (3.61): 


rx 2 pm g2 y 


I XI Jy 1 
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m ( g y 

dx 


'y i 
rv 2 


'Vl 


(x 2 ,y) 


dv 

dx 
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( w (x2,y) w (xi,y)) dy 


(3.62) 


where w is represents the angular velocity. If the patch is discretized along the edge, then 
we get the following approximation: 


rv 2 
'yi 


(^(x 2 ,y) [x i ■//) ) dy Ae ^ Wj 


(3.63) 


1+3 


where Ae is the length of the patch element and the summation is along edges 1 and 3 of 
the patch. Thus, from (3.62) and (3.63), we get 


fX 2 P y 2 02 

/ / d^ dxdy = Ae J2 Wi 

Jx 1 1+3 


Similarly, we can show 


px 2 py 2 Q2 y 
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Using (3.64) and (3.65) in (3.61) gives 
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(3.66) 


where ^ pzt denotes summation along all the edges of the patch. Equation (3.66) can be 
expressed in terms of complex amplitudes as 


T . , 24+nA a 2 V d 31 Y u (H + h)Ae^ j;(/ 
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(3.67) 
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where I and V denote the complex amplitudes of the PZT current and voltage respectively 
while Wi denotes the complex amplitude of the angular velocity at the i th element along the 
patch edge. The real part of the electric power is given by 

P e = ^Re(V*I) (3.68) 

where * denotes the complex conjugate. Using (3.67) and (3.69), we get 


Pe 






2 jjiM 

l-p p ) 


a 2 V*V d 31 Y n (H + h ) Ae - * ^ 

h 1 — u v 2 1 


(3.69) 


Note that for a thin patch, the moment per unit length (complex amplitude) can be expressed 
as a function of the applied voltage (complex amplitude) as 


M= d3lVu ^ + /j ) y 
1- u p 2 

The net moment applied to the patch element is, therefore, 


M 


A eM 


d n Y n (H + ft) 

1 — u p 2 


(3.70) 


(3.71) 


Using (3.71) in (3.69), we get 


P e = -Re 
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e 33 - 
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(3.72) 


Let 

£33 = ^ 33(1 ~ i&) (3.73) 

where £33 is the real part of the complex dielectric constant and S is the dielectric loss factor. 
Using (3.73) in (3.72), we get which simplifies to 


Pp = V| 2 e 33 d + ^Re 

pzt 


2 h 


(3.74) 
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Recall that the net input power (mechanical) due to the patch is given by 


P=^RelM*J2Wi) (3.75) 

\ pzt / 

Using (3.75) in (3.74), we get 

ft = vr^ + ft 

ft = ft + ft (3.76) 

where p£ is the dielectric loss. Equation (5.51) gives the relation between the net input 
electrical power and net input mechanical power. 


3.4.2 Real-time control algorithm to implement MAP for a single 
secondary actuator 

From Eq. (5.51), the real power in the secondary actuator can be expressed as 

Pse = Psi + Ps (3.77) 

where 

Pu= (3.78) 

is the dielectric loss, P se is the net electrical power in the secondary actuator, P s is the net 
mechanical input power from the secondary actuator and V s is the complex amplitude of the 
voltage applied to the secondary actuator. Since P s = 0 at optimal control, equation (3.77) 
shows that at optimal value of control, the net electrical power in the secondary actuator is 
given by the following baseline value: 


Pse = Pst (3.79) 

Equation (4.31) shows that setting P se — P S £ = 0 (with V s non-zero) provides a novel way to 
implement MAP control. 
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Next, we present simulations to illustrate the implementation of MAP control. The primary 
and secondary actuators are symmetrical with respect to the center of the plate. Specifically, 
(x p ,y p ) = (0.305,0.275) and (x s ,y s ) = (0.075,0.075). The voltage applied to the primary 
actuator is V p =10 volts. Figure 3.19 shows the real electrical power in the primary actuator, 
Ppe, with and without MAP control as a function of frequency. In case of perfect control, that 
is, P p = 0, the real electric power in the primary actuator is simply the dielectric loss. The 
top plot in Fig. 3.20 shows the variation of the total vibrational energy with the secondary 
voltage at the resonant frequency ui n = 620 Hz while the bottom plot shows the variation of 
real electric power in the primary and secondary actuators and the secondary loss. Because 
of symmetry in actuator location, the optimal secondary voltage is V s = 10 V. It can be 
seen clearly that the secondary voltage is optimal at the intersection of the secondary real 
electric power and the secondary loss as given by Eq. 4.31. 



Figure 3.19: Real electric power in the primary actuator with and without MAP control. 
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Figure 3.20: (top) Vibrational energy as a function of the secondary voltage at u> n = 620 Hz. 
(bottom) Real electric power in the primary and secondary actuators and the secondary loss 
as a function of the secondary voltage. 

3.5 Experimental Demonstration of MAP 


The theoretical findings just described provide the groundwork and guide an experimental 
validation. A test setup that closely mimics simply supported boundary conditions was built. 
One primary and one secondary patch actuator was symmetrically bonded to the plate struc- 
ture. The symmetrical location of the actuators helps validation as the optimal secondary 
voltage is expected to be equal in magnitude to the primary voltage. The theoretical findings 
show that at optimal control, the net input power from the secondary actuator is identically 
zero. The net input power from the secondary actuator can be indirectly computed by 
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measuring the real electric power supplied to the secondary actuator and calculating the 
dielectric loss. If the dielectric loss is accurately characterized, then MAP control can be 
implemented by solely monitoring the real electric power supplied to the secondary actuator. 


3.5.1 Overview of experimental verification 

We have shown that for a structure excited by a single primary actuator and controlled by 
a single secondary actuator, MAP control for vibration minimization of a single harmonic 
excitation can be practically implemented by monitoring, regulating and matching the real 
electric power of the secondary actuator to a predetermined value. Recall that at optimal 
control, the net input mechanical power from the secondary actuator into the structure 
diminishes. The net input power is related to the real electric power by the following simple 
relation: 

Ps = Pse — Pst 

where P s is the net input power, P se is the real electric power and P s £ is the electric loss. 
Thus, at optimal control, the real electric power in the secondary actuator ideally reduces 
to its electric resistive loss, or the dielectric loss for a piezoelectric actuator. The process 
of MAP control then simply reduces to the measurement and regulation of electric power , 
which is far easier than that of its mechanical counterparts, the vibratory angular velocity 
and bending moment. It is important to note that the dielectric loss is a known function 
of the piezoelectric material properties, excitation frequency, and the electric field applied 
(voltage). The dissipative power of the actuator is a quadratic function of the driving 
voltage. As the secondary actuator voltage increases, both the real power consumption and 
the electrical loss increase, although at different rates. Therefore, the optimal control voltage 
V s could be determined by the interception of the real time power measurement with the 
electrical loss, as illustrated in Fig. 3.20. For a given primary voltage V p — 10 volts, it shows 
variation of the power as the secondary voltage is increased from 0 to 15 volts. The optimal 
secondary voltage for the minimum vibrational energy is determined by the intersection of 
P se and P s i . The experimental implementation of the MAP control for each single frequency 
consists of two steps: 1) determining the actuator electrical loss as the control baseline and 
2) varying the secondary voltage with respect to the primary voltage until the measurement 
of the real power matches a predetermined resistive power loss. The MAP control of multiple 
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frequencies is obtained by matching the total electric power of all the sinusoidal excitations 
to the sum of the electric losses for all the frequencies. 


3.5.2 Test instrumentation 



Figure 3.21: Instrumentation of MAP experimental verification. 


The experimental process for MAP verification consists of multiple channel frequency sig- 
nal generating, voltage and frequency sweep, complex power measurement, and vibration 
monitoring, as shown in Fig. 3.21. The two audio power amplifiers with line transformers 
are used to drive the primary and secondary actuators having a common electric grounding. 
The real electric power in the secondary actuator is calculated from the secondary voltage 
and the current sampled by a precision resistor on the drive circuit. The real power of each 
actuator is monitored and recorded as function of drive frequency, voltage, and phase inputs. 
The vibration level of the plate structure is monitored by an accelerometer mounted on the 
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plate at a location away from most of the node lines of the modes (frequencies) being an- 
alyzed. The secondary power is then compared against the baseline electric loss prediction 
while the vibration level is monitored. The measurement system is programmed in Lab- 
view. The power consumption of the piezoelectric actuators as function of frequency is also 
independently characterized with an HP 4194 electric impedance analyzer. 


3.5.3 Structural test apparatus 

The experiment test apparatus is made of a simply supported thin aluminum plate. The 
plate is sandwiched between two steel plate frames, as shown in Fig. 3.22. The aluminum 
plate is V-notched around its four clamped edges to emulate flexible hinges with minimal 
rotational stiffness. The two steel frames are bolted together to create highly rigid mass 
blocks providing a nearly clamped boundary condition to the aluminum plate edges beyond 
the notch, as shown in Fig. 3.23. Two thin piezoelectric patches are bonded to the plate as 
the primary and secondary actuators. 



Figure 3.22: Test apparatus of simply supported plate for map validation. 
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Figure 3.23: The aluminum plate test apparatus with a primary actuator and a secondary 
actuator. 

The piezoelectric patch actuators are made of PZT 5H material with nickel plating on both 
side as electrodes. They are directly bonded to the aluminum plate with low viscosity 
adhesive Loctite 495. A uniform pressure is applied to the actuators by vacuum bag during 
the adhesive curing to minimize the thickness of the adhesive layer and minimize the electric 
and dielectric losses. Since the adhesive is semi-conductive, the two actuators therefore use 
the aluminum plate as a common electrode to complete their drive circuit. The dimensions 
and material properties for the plate and PZT patch are provided in Table 4.2. 


3.5.4 Implementation of MAP control 

A critical step for MAP validation is to establish the control baseline, which is to be derived 
from the electric power loss of the secondary actuator. We have surveyed the actuator power 
consumption characteristics in the frequency range of interest by an electric impedance 
analyzer. The HP 4194 electric impedance analyzer measures complex impedance as a 
function of frequency for low voltages (2.8 Vp-p). Figure 3.24 shows a frequency sweep of 
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Table 3.2: Geometry and material properties for plate and PZT patch. 


Variable 

Value 

Lx 

0.38 m 

Ly 

0.35 m 

H 

4.2e-3 m 

E 

70e9 Pa 

V 

0.01 

P 

2700 kg/m 3 

P 

0.33 

a 

0.05 m 

h 

0.13e-3 m 

Y u 

6.6el0 Pa 

di\ 

-190e-12 m/V 

^33 

1.53e-8 F/m 

V p 

0.35 


the real electric power deduced from the complex impedance measurement of both actuators 
mounted on the plate. 


1000 t 



Figure 3.24: Real electric power of the primary (P rp ) and secondary (P rs ) actuators, their 
dielectric losses at 10 V excitation. 
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As a result of the electromechanical coupling of the piezoelectric actuators to the plate, the 
real electric power plots of both actuators P p and P s reveal not only the electric power loss of 
the actuators themselves but also the structural dynamics and power coupling variation with 
frequency. A few observations can be made by this survey. The real electric power, linearly 
proportional to the mechanical power of the actuator, maps the plate vibration amplitude 
well. The actuator power increases substantially when the plate is resonance and collapses 
to a small baseline level when off-resonance. The actuator power off-resonance shows a 
general linear dependence with frequency. Although the primary and secondary actuators 
are symmetrical relative to the center of the plate and possess an identical geometry, their 
power consumption differs noticeably at resonance frequencies. This implies variation in 
their control authority when one is used as a primary and the other is the secondary. 

Figure 3.24 also compares the electric loss of the actuators as potential control baselines. 
One of them (Loss_C) is predicted by an analytical model (Equation 5.51) assuming a con- 
stant piezoelectric material loss tangent and the other (Loss_M) predicted by curve fitting 
the power measurement P rs and P rp at their off-resonance frequencies, respectively. The 
significant difference between these two loss estimates as control baseline candidates leads us 
to believe that a more comprehensive actuator characterization is necessary in order to val- 
idate the MAP control algorithm by the power interception approach mentioned early. The 
red dots indicate resonant frequencies and purple dots indicate the off-resonant frequencies 
selected for MAP validation later. The results are discussed in sections to follow. 


3.5.5 Actuator electric loss characterization 

We considered three methods to estimate the MAP control baseline from the actuator electric 
loss. 


1. Analytical linear relation in Equation 5.51 taking piezoelectric material loss tangent 
as constant (orange dash line in Fig. 3.24. This model, although simple and easy to 
follow, lacks sufficient fidelity when used in the power interception algorithm for MAP 
control. Measurement of the piezoelectric material by the electric impedance analyzer 
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reveals a variation of 250% in the loss tangent of PZT 5H piezoelectric in frequencies 
between 100 Hz and 2000 Hz, in constract to the nominal constants of 0.02 published by 
the material manufactures and academia as shown in Fig. 3.25. The loss tangent also 
largely depends on the electric field even at the low field intensity (0.2~2V/mil) and 
varies by more than 50% when the voltage changes by one order of magnitude (1~10V). 


2. Empirically curve fitting to the power measurement of the on-line actuator (mounted 
on the plate) at its off-resonance frequencies (orange solid line in Fig. 3.24). The base- 
line thus obtained is a direct measurement and a system calibration of the actuator on 
the structure represents the true power characteristics at all frequencies. 

3. Empirically curve fitting to an off-line or null actuator. The off-line actuator has the 
same size as the on-line but is mounted to a solid and thick aluminum block with very 
high flexural stiffness. When actuated, it induces nearly zero strain in the block so the 
patch actuator essentially behaves like a pure electric capacitor with only dielectric loss. 
Therefore, curve fitting to its power measurement represents a good control baseline 
of the optimal control condition, that is, when the vibration completely diminishes. 


The observed strong dependence on the material property and mismatch of the power in the 
initial MAP testing suggests that we may not be able to use the simple dielectric loss model 
for the control power estimate and would need to map out experimentally the actuator’s 
power loss in the entire frequency range of interest as a function of the applied electric 
field. We have constructed a null actuator, characterized its loss as a function of frequency 
and voltage, for comparison to the actuator on the plate. Figures 3.26 and 3.27 show the 
electric power survey of the plate-mounted actuator and null actuator, respectively. The 
power is measured between resonance frequencies of the plate, from 400 Hz to 2000 Hz, for 
drive voltages of 1V-10V. A large difference can be seen between the loss estimates for the 
two actuators. The null actuator consumes less power since its measurement was taken at 
nearly zero residual vibration and will likely represent the lower boundary of MAP control 
baseline. Comparatively, the plate-mounted actuator consumes more power for the same 
frequency and drive voltage, and therefore constitutes a poor estimate for control baseline 
as the vibration at the optimal control never truly becomes zero. 
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3.5.6 MAP control validation 


The MAP control of plate vibrations is performed using the instrument setup described 
previously in Fig. 3.21. The MAP validation is done by sweeping the secondary voltage 
while monitoring the plate vibration at various frequencies, both resonant and non-resonant, 
between 400 Hz and 2000 Hz. The secondary voltage at which the measured P s equals the 
control baseline P s -opt is then considered as the optimal voltage. For this study, the power 
loss measured with the actuator mounted on the aluminum plate was used. For the MAP 
control testing to follow, the baseline loss values at each resonance are calculated by linear 
interpolation between loss values computed at two adjacent non-resonance frequencies. The 
MAP control criterion was tested for validation on both single sinusoidal excitation and 
dual sinusoidal excitation. For most of the frequencies selected for testing, the MAP control 
criterion holds whether they are on-resonance or off-resonance. The degree of control author- 
ity demonstrated very much depends on the placement of the control actuator (secondary) 
relative to the source (or primary actuator). 



■■ ■■ Pre ■ Pim - loss tan Poly, (loss tan) 


Figure 3.25: Material dielectric loss of a null actuator mounted on a solid aluminum block. 
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Power of off-line actuator (null) @ off-resonance frequencies 


18 - 



Figure 3.26: The power loss of null actuator mounted on aluminum block as MAP control 
baseline metric. 


Power of on-line secondary actuator (mounted on plate) @ off- 
6 resonance frequencies 





Figure 3.27: The power loss of the on-line actuator mounted on aluminum plate as MAP 
control baseline metric. 

Single frequency: on-resonance and off-resonance 


Figure 3.28 presents the results of MAP control at four of the plate’s resonance frequencies; 
539, 816, 1136 and 1514 Hz, respectively, as illustrated by the red dots in Fig. 3.24. Each 


plot shows the primary, secondary power, the control baseline and vibration in dB as func- 
tion of secondary voltage V s while the primary voltage V p is fixed to 10 V. The interception 
of the P s with the control baseline is defined as the optimal control voltage V s -opt predicted 
by MAP. The secondary voltage that minimizes the plate vibration is V s -min dB. The dif- 
ference between the two voltages represents a deviation of the MAP prediction from the 
measurement. The case of 1136 ffz shows nearly perfect agreement between the prediction 
and the measurement. The primary and secondary power also meet at the optimal secondary 
voltage, reflecting the fact that the boundary conditions and actuator placement are nearly 
symmetric. The results at 539 Hz illustrate a good match between optimal and minimum 
plate vibration. However, the MAP result under-predicts the optimal control voltage at both 
812 Hz and 1514 Hz. While the vibration minimum occurs at 10 V and 9 V respectively, the 
secondary power intercepts the baseline at 8.6 V and 8.2 V, respectively. 

Figure 3.29 presents the results of MAP control at three frequencies between the plate 
resonance frequencies: 580, 1050 and 1450 Hz, as illustrated by the purple dots in Fig. 3.24. 
Each plot shows the primary, secondary power, the control baseline and vibration in dB as 
function of secondary voltage V s while the primary voltage V p is fixed to 10 V. As can be 
seen, the MAP optimal control voltage is a fair match to the voltage at minimum vibration 
for 589 and 1050 Hz. The prediction off-resonance overshoots the target, which is in contrast 
to the case on-resonance atl450 Hz, the MAP criterion fails as the secondary power does 
not intercept the baseline at all. 

The performance of MAP control is summarized in Fig. 3.30 which compares the measured 
control voltage (intersection of real electric power with the control baseline) to the measured 
minimum vibration voltage. As observed, the control error varies with frequency up to 15% of 
the targeted values. The MAP prediction at resonance frequencies appears to under-estimate 
the minimum vibration voltage, while the opposite happens between resonances. The overall 
deviation between the MAP optimal voltage and the voltage at minimum vibration averages 
at 4%. 
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Figure 3.28: Actuator power at four on-resonance frequencies of the plate. 


Multiple frequency excitation 


To validate the effectiveness of MAP control on a multiple frequency excitation, we have 
simultaneously applied two frequencies, 539 Hz and 812 Hz to the primary and secondary 
actuators. The analytical model of the multiple frequency excitation case shows that the total 
power and total loss (and therefore the control baseline) are the sums of that of each frequency 
components when the two frequencies are commensurate. Figure 3.31 shows summation of 
the primary power, secondary power, the control baselines by each frequency component. 


90 


Vibration. dB Vibration. dB 






(c) 

Figure 3.29: Actuator power at three off-resonance frequencies of the plate. 


Two resonances of the plate are chosen for this verification. The vibration reduction in this 
case is calculated by the change in total energy at both frequencies. The primary voltage is 
fixed to 5V while the secondary is swept from 1 V to 5 V for both frequencies. The MAP 
optimal control voltage is V s -opt= 4.2 V for both frequencies while the vibration minimum 
occurs at 4.5 V, as shown by Fig. 3.31. A total vibration reduction of 15 dB is observed at 
MAP optimal voltage. Consistent with the trend seen earlier at individual frequencies, the 
MAP voltage is below the minimum vibration voltage by 4.5%. 
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Figure 3.30: MAP prediction error analysis at on-resonance and off-resonance frequencies. 
V s -opt: predicted by MAP (interception), V s -min dB: actual corresponding to the minimum 
vibration. 
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Figure 3.31: The power and loss of primary and secondary actuators under two resonance 
frequency excitation, 539 Hz and 812 Hz, simultaneously. 


Table 3.3 summarizes the control effectiveness of MAP control at resonance and off-resonance 
frequencies. 
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Table 3.3: MAP control effectiveness at resonance, off-resonance and multi-frequency exci- 
tations. 


Mode 

On-Resonance 

Off-Resonance 

Dual Frequency 

Freq, Hz 

539 

816 

1136 

1514 

580 

1050 

1450 

539 + 816 

MAP Prediction Error % 

-6 

-15 

-2 

-8 

6 

7 

n/a 

-4.5 

dB in reduction 

-33.7 

-47.2 

-39.1 

-42.2 

-1.7 

-7 

n/a 

-20 


The experimental investigation of the simply supported plate with a single primary and 
single secondary actuator has demonstrated the effectiveness of the proposed MAP control 
algorithm for active structural vibration control. It has also identified some of the potential 
challenges with the implementation. Using an empirically derived actuator electric loss 
model as the control baseline, the MAP control approach produces good vibration reduction 
in a broad frequency range of the structure except at one tested frequency. The difference 
between accuracy of the MAP optimal voltage and the voltage corresponding to minimum 
vibration varies by frequency and modal status (resonance or non-resonance) of the structure 
with observed average difference below 10%. The difference between the MAP voltage and 
the minimum vibration voltage depended on whether the excitation frequency was between 
resonances or at a resonance. 

Vibration reduction produced by the secondary actuator is very effective at structural reso- 
nance mode for both single frequency and dual frequency vibration (> 20dB). However, it is 
much less effective between resonances off-resonances (> 1.7 dB). This does not necessarily 
suggest that MAP control is not effective at off-resonance frequencies. The effectiveness of 
MAP control depends on the location of the secondary actuator with respect to the primary 
actuator. Optimizing the location of the secondary actuator can substantially improve the 
performance of MAP control between resonances (see §3.3.6). 

The most critical step in experimental implementation of the MAP is identified as the esti- 
mate of the control baseline with the specific piezoelectric actuator on the structure. It is 
essential that the electric loss metric be based on or calibrated against an actual measurement 
of the actuator being used for control. 
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3.6 Summary 


We summarize this chapter by stating the key accomplishments 


• The theory of MAP was developed for a structure excited by a single primary excitation 
and controlled by a single secondary excitation. 

• An analytical model for a simply supported plate excited and controlled by piezoelectric 
actuator was developed. 

• An electromechanical coupling model of a piezoelectric patch actuator on a structure 
is developed and integrated with the MAP model. 

• MAP AVC was demonstrated on a simply supported plate experimentally at discrete 
excitation frequencies, both resonant and off- resonant. 
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Chapter 4 


MAP Control for Multiple Primary 
and Multiple Secondary 


In chapter §3, we laid the foundations of MAP control and developed the theory for single 
primary and single secondary actuator. Simulations and experimental results were presented 
to validate the theory. In this chapter, we focus on extending the theory of MAP control to 
multiple primary excitations and multiple secondary actuators with an emphasis on devel- 
oping modeling capability and a real time control algorithm to demonstrate the technology’s 
effectiveness and ease of implementation on a rotorcraft interior panel structure. In §4.1, we 
develop the system model for MAP control of a simply supported plate that is excited by 
multiple primary piezo-electric (PZT) patch actuators and controlled by multiple secondary 
PZT patch actuators. In §4.2, we present numerical simulations for a simply supported 
plate excited by two primary patches and controlled by two secondary patches and evaluate 
MAP control performance as a function of the number and location of secondary patch ac- 
tuators. In §4.3, we report experimental results that demonstrate MAP control of multiple 
primary excitations with multiple secondary actuators and quantify MAP performance and 
effectiveness in terms of vibration reduction and secondary power requirements. 
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Figure 4.1: (a) Schematic of a simply supported rectangular plate excited by PZT actuators 
(b) Distributed moment generated by a patch. 

4.1 System Model For Theoretical Simulations 


4.1.1 System description 

We consider a simply supported rectangular plate that is excited by N p primary PZT patch 
actuators and controlled by N s secondary PZT patch actuators. For sake of illustration, Fig 
4.1(a) shows the plate excited by one primary PZT patch and one secondary PZT patch. 
Each actuator is driven by a voltage source. The primary patch acts as the disturbance while 
the secondary patch is the control actuator. Let L x , L y and H be the length, width and 
thickness of the plate respectively. The PZT patches are bonded to the plate surface and are 
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assumed to be square with width a and thickness h. When an electric field is applied to the 
patch in the z direction, the patch expands and contracts in both the x and y directions. The 
resulting strain field generates a distributed moment along the edges of the patch, labeled 1 
to 4, as shown in Fig 4.1(b). 

We make the following assumptions: 


• The plate is isotropic. 

• The patch is thin and the patch inertia and stiffness are ignored. 

• The primary and secondary patches are identical. 

• No losses are added due to the patch bond to plate 

• The boundary conditions are simply supported. 


4.1.2 Total input power minimization for multiple primary and 
secondary patch excitation 

In the Multi- Actuator Data package submitted previously this year, we developed the theory 
of MAP for multiple primary sources controlled by multiple secondary sources. For simplicity, 
we assumed the primary and secondary sources to be point forces. However, the experimental 
setup uses PZT patches as actuators that generate line moments along the patch edges 
instead of a point force. For completeness, we present the theory of MAP for multiple 
primary and multiple secondary patch actuators to match the experiments. It should be 
noted that the analysis and some key observations are similar to the point forces previously 
reported but it is important to close the loop on the theory and experiments. The point 
forces are replaced by moments and the linear mobilities replaced by equivalent angular 
mobilities. 

The theoretical development of input power minimization for multiple primary and multiple 
secondary actuators follows closely the analysis for a single primary actuator controlled by 
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single secondary actuator. Let M pi and M si be the moments per unit length (complex 
amplitude) generated by the i th primary and i th secondary patch respectively. Let the 
primary and secondary patch edges be discretized into an equal number of elements and let 
Ae be the length of the edge element. Then, the moments generated by each element of the 
jth p r j mai y anc [ 2 th secondary patch, respectively, are 

M pi = Ae M pi 

M si = A eM si (4.1) 

Let 

M p = [M p i M P 2 . . . M p atJ (4.2) 

represent the moment vector due to N p primary patches and 

Mj = [M s i M s2 ... M sNs ] (4.3) 

represent the moment vector due to N s secondary patches. We assume that the primary and 
the secondary patch moments are harmonic with frequency u>. Let 

M t = [Mj Mj] (4.4) 

be the total moment vector. The total power input P is simply the sum of the power inputs 

from all the primary and secondary patches. That is, 

P = P p + Ps (4.5) 

where P p and P s are the total power inputs from all the primary and secondary patches 
respectively. The total input power can be expressed as 

P = e{M H PM} (4.6) 

where P is the equivalent rotational mobility matrix. The reader is referred to §B for a 
detailed derivation of (4.6) and the analytical expression for P. P is symmetric due to 
reciprocity. In this case 

P = ^M H Re{P}M = ^ M h SM (4.7) 
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where S is the real part of P. Decomposing into primary and secondary patch components, 
we get 

M t = {M p T Mj} (4.8) 

where M p and M s are the complex amplitudes of the primary and secondary patches re- 
spectively. Similarly, S can be decomposed as 

(4.9) 

where S sp = S ps . Using (4.7), (4.8) and (4.9), we can express the total input power as 

P = ^ [Mp SppM p + S ps M s + M^S sp Mp + M?S ss M s ] (4.10) 

Note that P is a quadratic function of M s . Minimizing P with respect to M s , we get the 
optimal secondary patch moments are [41] 

m; = -S-'S, P M P (4.11) 

The optimal secondary patch moments depend not only on the magnitude and phase of the 
primary patch moments but also on location of the primary patches with respect to the 
secondary patches. In general, the primary patch moments are unknown. In what follows, 
we assume that we know the location of the primary patches. However, we do not know the 
magnitude and phase of the primary patch moments. 



The net input power from the secondary patches is 

P s = e{M»P sp M p + MfP ss M s } (4.12) 

Since P ss is symmetric, we get 

R e(M?P ss M s ) = M^Re. ( P ss ) M s = Mf S ss M s (4.13) 

Using (4.13) in (4.12), we get 

Ps = e{M*P sp Mp} + l -M?S ss M s (4.14) 

We consider the following two cases. 
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Primary patch moments in-phase or 180° out-of-phase 


If the primary moments are in phase (or 180° out-of-phase), then 

M p = M p e i<p => M f = M^e~^ (4.15) 

The real input power from the secondary patches at optimal condition is 

P; = 1r e{-M r H S r ,S„ 1 P, r M r } + ±M»S r ,S;, 1 S„S 7 , 1 S, r M r 

= -jR e{MlS p ,S^P, r M r ] + 

= -- MlS v ,S-}S,,M v + \M T r S r ,S-}S, v M r 

= 0 (4.16) 

Thus, the real secondary power is zero when the the primary moments are either in-phase 
or 180° out-of-phase. 


Primary patch moments with arbitrary phases 

In this case, the net input power from the secondary patches at optimal condition is 

P’s = -lR e{M«S p ,S- 1 P, p M p }+ 1 -MfS p ,S;, 1 S, p M p (4.17) 

Note that P* = 0 when S ps S~g P sp is symmetric. If 5 , ps 5 , “ 1 P sp is symmetric, then we have 
Relief SpsS-'PspMp} = M p H Re{Sp S S^P S p}Mp 

which gives 

P’s = -\MfRe{Sp,S- 1 Ps P }M r + 1 -MfS p ,S;, 1 S, p M p (4.18) 

= -\M»S p ,S^S, r M p + 1 -M»S r ,S-, 1 S, p M r (4.19) 

= 0 (4.20) 
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(4.21) 


Now, for S ps S s J-P sp to be symmetric we need 


c <j-i p _ p e-i c 1 

^ps^ss ± sp J ps^ss 


Using P sp = S sp + jX sp , we get 


Sp I J Spt 

which gives 

Equation (4.23) holds if 


S ps S-\S sp +jX sp ) = ( S ps ■ jX ps )S~ 1 S 


sp 


<? _1 V — V <?-! <? 

^ps^ss y *-sp — ^ps&ss °sp 


Ssp (3X sp 


(4.22) 

(4.23) 

(4.24) 


In order to better understand the above result, consider the case where there are two primary 
patches and one secondary patch. In this case, 


Psp [‘S'spi 4“ ]X Spi S S p 2 jXgp^l 


We get 


c e-i p _ 

k -’ps J ss Jr sp — Q 

dee 


s. 


[‘S'spi 4” j-^-Spi S sp2 “4 jX sp2 ] 


‘-’sp 1 

q 

u sp2 

S spi 4~ j S spi X spi S spi S sp2 j S spi X S p 2 

S S p 1 S sp2 + jS sp2 X spi S sp2 + jS sp2 X sp2 


From (4.26), for S ps S ss l M sp to be symmetric, we get 


S, P1 X,„ = S, P2 X, P , => = = ,8(3 ay ) 


X S p , X S p 2 


Equivalently, 


Ssp [ —‘sp I Sgp 2 ] 

= [ (3X SP1 (3X SP2 ] 

= px sp 


(4.25) 


(4.26) 


(4.27) 


(4.28) 
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Equation (4.28) puts a constraint on the location of the secondary patch with respect to the 
primary patches in order to have the real secondary power equal to zero. It implies that 

ipu = -012 (4.29) 

where 0n = arg(P spi ) and 0i2 = arg(P SJ)2 ). That is, the phase introduced due to the 
equivalent rotational cross-mobility is the same for both the primary-secondary pairs, ( p\ — s) 
and (p 2 — s). 


4.1.3 Net electrical power in secondary actuator at optimal con- 
trol 

In an actual experiment, the net input (mechanical) power from the secondary patch actua- 
tors is hard to measure. Rather, one can easily measure the electrical power by monitoring 
the voltage and current in the patches. The relation between the net electrical power and 
the net input power for the i th secondary is given 

P e „ = \ i ^-\ v 4'^S + P,i (4.30) 

where P ^ is the net electrical power in the i th secondary patch, V si is the complex amplitude 
of the voltage applied to the secondary actuator, e 33 is the real part of the complex dielectric 
constant and 6 is the loss factor. See §3.4 for detailed derivation of the relation between 
net input power and net electrical power for a generic PZT patch. Since P S j = 0 at optimal 
control, equation (4.30) shows that at optimal control, the net electrical power is equal to 
the dielectric loss in the patch: 


PI = - 


lu a 2 


2 h 


|fdsi| 2 ^33^ 


(4.31) 
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4.2 Numerical Simulations 


Here, we present simulation results for a simply supported plate that is excited by two 
primary patch actuators and controlled by one or more secondary patch actuators. The 
effectiveness of MAP control depends not only on the position of the secondary patches with 
respect to the primary patch actuators but also on the number of secondary patch actuators. 


4.2.1 Simulation parameters 


y 
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Jbpi 


primary 


Vpi 



X x 
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Figure 4.2: Schematic of the plate system for simulation. 

A schematic of a simply supported plate system for the simulation studies is shown in Fig. 
4.2. For sake of illustration here, one primary and one secondary patch actuators are shown. 
The coordinates of the i th primary and secondary patch centers are ( x pi ,y pi ) and ( x si ,y si ) 
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respectively. We assume that the plate is made of aluminum. Table 4.1 lists geometrical and 
material properties for the plate and PZT patches. 

Table 4.1: Simulation Parameters. 


Variable 

Value 

Lx 

0.38 m 

Ly 

0.35 m 

H 

4.2e-3 nr 

E 

70e9 Pa 

V 

0.01 

P 

2700 kg/ nr 3 

P 

0.33 

a 

0.05 nr 

h 

0.13e-3 nr 

Y n 

6.6el0 Pa 

d'M 

-190e-12 nr/V 

^33 

1.53e-8 F/nr 

U p 

0.35 

5 

0.022 


4.2.2 MAP Performance as a Function of Location of Secondary 
Actuators 

The primary focus of this study is to evaluate MAP performance as a function of the num- 
ber and location of secondary actuators for a given set of primary actuators. We assume 
that the primary patch voltages are either in phase or 180 degrees out of phase. We first 
consider the plate being excited by two primary patches and controlled by a single sec- 
ondary patch. Figure 4.3 shows the total power input with and without MAP control for a 
particular configuration of the primary and secondary actuators. At resonance frequencies, 
the total power input is proportional to the global vibrational energy. Therefore, minimiz- 
ing the total power input is equivalent to minimizing the global vibrational energy. The 
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locations of the primary and secondary actuators are arbitrarily chosen. In this example, 
(x p i,y p i) = (0. 51^/4, 0. 51^/4), ( x p2 ,y p2 ) = (3.5L x /4, 3.5L y /4) and (x sl ,y sl ) = (L x /2,L y /2). 
The primary patch voltages are V p i = 10 V and V P 2 = 10 V. It can be seen that while some of 
the resonance peaks are attenuated, some are not. This depends on the position of the sec- 
ondary actuator with respect to the primary actuators. Consider two more cases as shown in 
Figures 4.4 and 4.5. In Fig. 4.4, (x pl ,y pl ) = {2L X /A, 0.5L y /4), (x p2 ,y p2 ) = (2L x /4, 3.5L y /4) 
and (x a i,y 8l ) = (L x /2,L y /2) while in Fig. 4.5, (x pl ,y pl ) = (0.5L x /4, 0.5L y /4), ( x p2: y p2 ) = 
(3.5L x /4, 0.5L y /4) and (x s i,7/ s i) = (L x /2, L y /2). In both cases, V p \ = 10 V and V p2 = 10 V. 
From Fig. 4.4, it can be seen that given the location of the primary actuators, the secondary 
actuator location is optimal for suppressing all the resonance peaks in the 500-2000 Hz fre- 
quency range. However, for the primary actuator locations shown in Fig. 4.5, the location 
of the secondary actuator is not optimal. Some of the resonance frequencies are attenuated 
while some are uncontrolled. In general, one may need more than one secondary actuator 
for improved MAP performance over a broad range of frequencies. 
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Figure 4.3: Plate excited by two primary (red) actuators and controlled by a single secondary 
(blue) actuator. Total power input with and without MAP control shown on right. 
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Figure 4.4: Effect of primary actuator locations on MAP control. The secondary actuator 
location is optimal for suppressing resonance peaks in the 500-2000 Hz frequency range 


Primary (red) Secondary (blue) 




Figure 4.5: Effect of primary actuator locations on MAP control. As the primary actuator 
locations are varied, the given secondary location is not optimal for suppressing all the 
resonance peaks in the 500-2000 Hz range. 

4.2.3 MAP Performance as a Function of Number of Secondary 
Actuators 


It is expected that MAP performance will improve as the number of secondary actuators is 
increased. We fix the location of the primary actuators as shown in Fig. 4.5. We start with 
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Figure 4.6: MAP control with two secondary patch actuator. 


Primary (rad) Secondary (blue) 




Figure 4.7: MAP control with three secondary patch actuators. Adding one additional sec- 
ondary actuator improved MAP performance at a given frequency while reducing vibration 
at other frequencies too. 

one secondary actuator and add one extra secondary actuator each time. Figures 4.5-4. 7 
show the configuration with one, two and three secondary actuators respectively. It can 
be clearly seen from figures 4.5-4. 7 that MAP performance significantly improves as more 
secondary actuators are added. The resonance frequencies that were not controlled by one 
secondary actuator, as shown in Fig. 4.5, are effectively controlled by adding more secondary 
actuators. 
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Frequency (Hz) 


Figure 4.8: MAP control with three secondary patch actuators with a different secondary 
actuator configuration. 


Given a fixed number of secondary actuators, MAP performance will depend on the location 
of the secondary actuators. Optimal location of the secondary actuators may also depend 
upon the location of the primary actuators and the frequencies that need to be controlled. 
This can be seen from Eq. 4.11 that shows the dependence of the optimal secondary moment 
on the cross angular mobility S sp which depends on the location of the secondary actuators 
with respect to the primary actuators. This is illustrated by comparing the MAP perfor- 
mance for two different secondary actuator configurations, as shown in Figures 4.7 and 4.8. 
One of the secondary actuators is moved to a different location. It can be seen that changing 
the location of one of the secondary actuators in blue made some of the resonance peaks, 
otherwise controlled, uncontrollable and deteriorated MAP performance. The location of 
the secondary actuators can be optimized by maximizing the sum of the mode shapes of the 
frequencies to be controlled. 


4.2.4 Secondary Input Power at Optimal Control 

When the primary actuators are in phase or 180 degrees out of phase, the total secondary 
input power is zero at optimal control regardless of location of the secondary actuators. If 
the primary actuators have arbitrary phase, then the secondary location has to satisfy a 
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constraint, as given by Eq. 4.29, for the secondary power to be zero at optimal control. 
Moreover, if there are more than one secondary actuators, then the input power from each 
secondary actuator is zero at optimal control. This important observation is key to devel- 
oping a stable iterative algorithm to find the optimal secondary voltages. Note that the 
stable iterative algorithm can be developed even when the primary actuators have arbitrary 
phasing. Figure 4.9 shows the input powers from the primary and secondary sources for the 
actuator configuration shown in Fig. 4.5 at two different excitation frequencies oj = 826 Hz 
and u j =724 Hz. It can be seen that the total secondary input power (in this case, there 
is only one secondary actuator) is identically zero when the total input power is minimum. 
Note that a is a scaling factor defined as a = V sl j V* for all i where V si is the voltage applied 
to the i th secondary actuator and V* t is the optimal value of the i th actuator. Thus, a = 0 
implies no control whereas a = 1 implies optimal control. In this case, V p i = V P 2 = 10 V. 
V* x = 7.058 V and V* x = 7.06 for oj = 826 Hz and oj = 724 Hz respectively. 

Figure 4.10 shows the input powers from the primary and secondary sources for the actuator 
configuration shown in Fig. 4.6 at two different frequencies: u j = 800 Hz (off-resonance) 
and oj =724 Hz (resonance). Since there are two secondary actuators, we plot the input 
powers from either actuator along with the total input power and the primary input power. 
It can be seen that the secondary input power from either actuators is identically zero when 
the total input power is minimum. In this case, V p i = 10 V and V P 2 = 5 V. (17*, , V* 2 ) = 
(5.4395, 0.7616) V and (V s \ , V s * 2 ) = (4.0181, 1.8153) V for w = 800 Hz and u = 724 Hz 
respectively. 


4.2.5 Effect of High Damping on MAP Control 

For simplicity, we consider one primary excitation controlled by one secondary patch ac- 
tuator. We first consider low structural damping with r] = 0.01. Figure 4.11 shows the 
secondary electrical power, the secondary loss and the total input power as a function of 
the secondary voltage. The primary voltage is = 10 volts and the excitation frequency is 
oj = 620 Hz. It can be seen that for this case of low damping, the intersection between the 
secondary electrical power and loss curves is clear and easy to detect. 
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[Vpl Vp2]=[10V 10V] (0=826 Hz 


[Vpl Vp2]=[10V 10V] (u=724 Hz 




(a) (b) 

Figure 4.9: Total, Primary and Secondary Input Powers for the configuration shown in Fig. 
4.5. (a) c o = 826 Hz (b) cj= 724 Hz. 
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Figure 4.10: Total, Primary and Secondary Input Powers for the configuration shown in Fig. 
4.6. (a) lu = 800 Hz (b) u = 724 Hz. 
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Figure 4.11: Intersection of secondary electrical power and loss with damping equal to 1%. 
The primary voltage is V p = 10 volts and the excitation frequency c o = 620 Hz. 

As damping is increased without increasing the primary excitation level, the secondary elec- 
trical power approaches the loss curve. This makes determination of the intersection point 
difficult. This is illustrated in Fig. 4.12(a) where the damping is artificially increased to 
8 %. The primary excitation level was kept at V p — 10 volts. In this case, the primary 
source is not strong enough to induce the excitation level sensed by the secondary actuator 
for implementing MAP control. However, if the primary excitation level is increased, enough 
excitation will be sensed by the secondary actuator. The separation between the secondary 
electrical and loss curves will increase resulting in a clean intersection. This is illustrated in 
Fig. 4.12(b) where the primary excitation voltage is increased from 10 volts to 30 volts for 
the same high damping value of 8%. 


Ill 




(a) (b) 

Figure 4.12: Intersection of secondary electrical power and loss with damping equal to 8%. 
(a) V p — 10 volts and (b) Increasing values of V p . 

4.3 Experimental Demonstration of MAP Control 


This section describes experimental validation of MAP control for multiple primary and 
secondary actuators. In §4.3.1, we first describe MAP validation using impedance head 
measurements to demonstrate how the secondary input power is zero when the total input 
power is minimized. Impedance head measurements allow us to calculate input powers 
directly. In §4.3.2, the testbed for real-time MAP implementation is described. Finally 
in §4.3.3, experimental results for MAP control with two primary and two secondary PZT 
actuators are presented. 


4.3.1 MAP verification with Impedance Head Measurements 

Minimum actuation power is verified by measuring the actual mechanical powers that are 
injected into the vibrating structure by primary and secondary sources. This verification 
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step demonstrates the foundations of MAP without considering piezoelectric actuators or 
how the mechanical power in the secondary sources is sensed. The mechanical power injected 
into a structure is easily computed using an impedance head that measures simultaneously 
acceleration (or velocity) and force according to 


The impedance head at both primary and secondary sources allows for the estimation of the 
total mechanical power injected into the structure by simple summation of the individual 
mechanical powers at each source. This assumes most if not all the mechanical power is 
contained in flexural waves as measured by the impedance heads. 


Impedance Head Experimental Setup 

Verification of MAP using impedance head was performed on a freely suspended aluminum 
plate. The plates dimensions were 15.1 inches x 14.9 inches x 0.125 inches. Impedance 
heads (PCB model 288D01) were mounted onto the structure at both primary and sec- 
ondary sources. A picture of the test setup is presented in Fig. 4.13. Free-floating shakers 
injected mechanical power through the impedance heads, and normal to the plate surface. 
A Lab VIEW code measured acceleration and force at all impedance heads and calculated 
the respective mechanical power. Single frequency excitations were generated and fed to the 
sources. With a set amplitude and phase of the primary source(s), the phase and amplitude 
of the secondary source were varied until MAP was achieved. Total mechanical power in the 
structure was calculated by simple summation of the individual powers of all sources. 



(4.32) 
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Figure 4.13: Plate Instrumented with Impedance Head (shakers not shown). 


Results 


The measured mechanical power for a single primary and single secondary sources when the 
secondary phase and amplitude are varied are presented in Fig. 4.14. The primary excitation 
was at 843Hz, a plate resonant frequency. Locations of the primary and secondary sources 
ensured both sources coupled well with this structural mode. The global structural power is 
also presented. Due to the aluminum plates very low damping, the phase of the secondary 
source is approximately either 0 or 180 degrees out of phase compared to the primary source. 
In the figure, phase information is contained in the sign of the amplitude voltage, with a 
negative value representing a 180 phase angle. For this particular source configuration, at 
MAP, the secondary has is 180 degrees out of phase compared to the primary. As expected, 
as gain is increased (becomes more negative in this case), the secondary’s real power becomes 
negative before increasing again, becomes 0 where the global power is minimized, and then 
becomes positive where the total power in the structure increases. This verifies analytical 
calculations as presented in §4.1.2. 
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Secondary Voltage [V] 

Figure 4.14: Impedance Head Measurements of Real Mechanical Powers. 


While impedance heads measured point force and acceleration, the theory also holds for 
moment excitations. This test setup provides a platform for the simultaneous development of 
MAP using piezoelectric actuators, a methodology that uses electrical quantities to estimate 
mechanical power without any impedance heads, and the development of control strategies 
with multiple primary and secondary sources and using MAP as a control objective. 


4.3.2 MAP Testbed and Real-Time Implementation 

We have developed control software and an instrument system that enables semiautomatic 
MAP implementation and allows for testing of various primary and secondary excitation and 
sensing. In order to implement MAP control to real structures without knowledge of the 
primary excitation and to facilitate more efficient tests with different configurations, it was 
considered important to establish a flexible hardware setup and to develop comprehensive 
software modules that automate the entire procedure. The control software and test system 
was upgraded to control four actuators simultaneously. Sections 4.3.2 and 4.3.2 contain de- 
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tailed descriptiosn of two important software modules, a control objective function estimator 
and a frequency detector. In §4.3.2, integration of the newly developed software modules 
with the main MAP control software is described in detail. 


Determination of Actuator Aggregated Loss 


A software module was created to estimate aggregated loss of the secondary actuator in- 
cluding dielectric loss of the PZT and loss due to the residual vibration of the structure left 
from the control. The software module is shown in Fig. 4.15. This module was developed in 
Lab VIEW and allows users to select voltage and frequency ranges before actuator excitation. 
The software has two loops, one for voltage sweep and the other for frequency sweep. First, 
the excitation voltage is fixed to a starting voltage, and the frequency is swept within the 
specified range. As frequency increases, the real power at each frequency is calculated di- 
rectly from voltage at a power amplifier and current at a resister connected to the actuator. 
After this step, voltage is increased and the frequency sweep is repeated. Finally, a raw 
power loss map is generated and stored as an array variable and sent to a two dimensional 
regression module for the aggregated loss estimation. 

The two dimensional regression module provides curve fitting to the power measurement of 
the actuator. The embedded algorithm searches for local minima in the measurement which 
are induced by off-resonance frequencies. Then, interpolation between the fitting points is 
performed to create a complete control objective function with respect to frequency and 
voltage. A heavy dependence of material properties and power mismatch were observed 
in the initial MAP testing. This stimulated the need for mapping out actuator power loss 
experimentally over the entire range of frequencies and electric field instead of using a simple 
dielectric loss model. The 2D map is stored as a spread sheet and read when the MAP control 
is performed. 
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Figure 4.15: Actuator electric loss characterization. 


Detection of Frequency and Structural Mobility Phase 


A software module was created to use the secondary actuator (control actuator) as a passive 
sensor to determine frequencies to suppress. In the experimental setup, a primary function 
generator creates a continuous sinusoidal waveform function with respect to f p , V p , and (ft p 
specified by a user. Then, the current through a resistor connected to the PZT is read by 
a digitizer. A built-in Lab VIEW function then computes the Fast Fourier Transform of the 
current signal to determine dominant frequency components in the signal. The FFT and 
peak frequency detection thus allow the primary frequency f p to be determined without 
using primary actuator information. The frequency of the secondary actuator, f s , is then 
set equal to f p . 

The cross mobility phase between the primary and secondary actuators is determined by 
comparing the phase between V p and I s . This is done in the off-line plate characterization 
prior to the MAP control. It can also be extracted from plate’s analytical models. The 


117 




Figure 4.16: Frequency and mobility phase detection for mMAP control. 


mobility phase detection technique is empirical and relies on experiments. Note that the 
cross-correlation function in Fig. 4.16 is also a built-in Lab VIEW function that provides the 
phase at a select frequency between two time domain signals. When piezoelectric transducers 
are used, the phase of I s can be changed by 7t/2 due to capacitance of the transducers. In 
this case, an actuator phase correction is needed. f s and mobility between the primary and 
the secondary sources are the output of this module to be used for the final MAP control. 


MAP control independent of the primary excitation 


The experimental process for MAP verification consists of multiple channel signal genera- 
tion, acquisition and vibration monitoring as shown in Fig. 4.17. The updates to the test 
setup include: 1) hardware update to enable multi-actuator excitation and sensing, and 2) 
the development of a software module that incorporates objective function acquisition and 
frequency detection to real-time MAP control. The overall measurement system has been 
developed based on Lab VIEW 8.6 software that controls an NI PXI chassis. The chassis con- 
tains analog to digital converters and waveform generators connected to four audio power 
amplifiers with line transformers. The real power of each actuator is monitored and recorded 
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Figure 4.17: Frequency and mobility phase detection for MAP control. 


as a function of drive frequency, voltage, and phase with respect to phase reference defined 
by the primary actuator. The real electric power in the secondary actuator is calculated from 
the secondary voltage and the current sampled by a precision resistor on the drive circuit. 
The vibration level of the plate structure is monitored simultaneously by an accelerometer 
mounted on the plate at a location away from most of the node lines of the modes being 
analyzed (see Fig. 4.18). 

Figure 4.17 illustrates how a secondary actuator is controlled to track the optimal control 
frequency, phase, and voltage. The objective function for the secondary actuator is the sum 
of electrical loss and mechanical loss (aggregated loss) as a function of frequency and voltage 
and is predetermined. A frequency detector analyzes signals measured by the secondary 
actuator and determines frequency contents of the disturbance in the structure f s . Then, the 
software determines optimal phase (</>,,) for the secondary actuator. The optimal secondary 
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voltage is determined at two different secondary phases, ip s and %jj s + 90 degrees. The ratio of 
the two voltage values is used to estimate (ft s . Once (f> s is determined, the secondary voltage 
is swept again to find the optimal control voltage. The developed software is designed to 
conduct frequency, phase and optimal control voltage estimation with a few clicks by users 
in pursuit of fully automated real-time MAP control. 


4.3.3 Experimental Validation of MAP Control with Multiple Ex- 
citations 

The principle of MAP for single primary (excitation) and single secondary (control) has been 
experimentally validated and demonstrated via piezoelectric collocated sensor and actuator 
on a simple metal plate. This year, we have extended MAP control theory to multiple 
primary excitations and multiple secondary actuators (mMAP) and validated it through 
the enhanced control software and hardware platform with emphasis on the phase control 
and optimization. We have also initiated MAP demonstration on a more realistic helicopter 
structure and with high damping. 

Theoretical investigations conducted early in Year 2 identified that: 1) principle of MAP is 
valid for mMAP and vibration suppression is equally effective if the primary and secondary 
actuation configuration can meet certain symmetry conditions, 2) MAP control implemen- 
tation algorithm through gain control (magnitude optimization) on single excitation is also 
applicable to the multiple excitation, and 3) for a more realistic excitation where arbitrary 
phasing of the primaries is dominant, the secondary phase optimization becomes an essential 
and necessary part of mMAP control. 


Description of mMAP validation 


The experimental setup described in $ 3.5 was used for the mMAP investigation and de- 
velop fundamental understanding. To accommodate the mMAP test, two piezoelectric patch 
actuators are added and the four actuators are now symmetrically located on four corners 
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Figure 4.18: The aluminum plate test apparatus with two primary and two secondary actu- 
ators. 

respectively to meet the nominal symmetry requirement of mMAP control. This placement 
also facilities flexibility in pairing actuators to form various symmetric and asymmetric con- 
figurations, see Fig. 4.18. Also depicted in this figure is the accelerometer used to sense 
local vibration. The two new actuators (P 2 an d S 2 ) are bonded to the aluminum plate with 
direct electrodes in back as the common grounding to minimize the impact of the adhesive 
between the actuators and the plate on the actuator aggregated power loss. This is critical 
to determining optimal control of the secondary actuators. The dimensions and material 
properties of the plate and PZT patch are provided in Table 4.2. 

The process of mMAP control (with gain only) is similar to the single MAP except that only 
one variable is swept at a time while the others are kept constant in magnitude and phase. 
The control objective is the actuator aggregated loss for the two control actuators. These 
losses are predetermined in the form of a 2D matrix of power loss as a function of frequency 
and voltage. 

In searching for the optimal control voltage of the secondary actuators, the two primary 
actuators excite the plate with a constant voltage V p at select frequencies. The voltage of 
the first secondary actuator, V s i, is swept at a constant frequency while comparing its power 
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Table 4.2: Geometry and material properties for plate and PZT patch. 

Variable Value 


Lx 

0.38 m 

Ly 

0.35 m 

H 

4.2e-3 m 

E 

70e9 Pa 

V 

0.01 

P 

2700 kg/m 3 

P 

0.33 

a 

0.05 m 

h 

0.13e-3 m 

Y u 

6.6el0 Pa 

d;n 

-190e-12 m/V 

^33 

1.53e-8 F/rn 

U p 

0.35 


to its actuator aggregated loss until P s i, the power in secondary actuator, equals the loss 
power, at which point vibration should be a minimum. The corresponding voltage is 
the estimated optimal MAP control voltage for the first control actuator. The next step 
is to excite this secondary actuator at a constant V* x while sweeping the voltage input to 
the second control actuator voltage V s2 until P s2 equals its aggregated loss at V S 2 =V* 2 . The 
vibration of the plate is monitored by an accelerometer in the center of the plate. Applying 
V s * and V* 2 to the control actuators Si and S 2 , respectively, further reduces the vibration 
to a new minimum. A few iterations between the two control actuators may be required to 
converge the control inputs to their global optimal values and further reduce the vibration. 


Multiple Primary Excitation with Simple Phases 


As the first example, we demonstrate performance of MAP with two primary excitations 
and two secondary controls on the lightly damped aluminum plate. Figure 4.19 shows the 
optimum search of the second control actuator with both primary actuators Pi and P 2 on 
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and the first control actuator Si excited at a fixed V* 2 . With both control actuators at the 
optimal control, the vibration is reduced by 20 dB. 
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Figure 4.19: Optimum search of mMAP at 1100 Hz mode. The purple vertical dashed lines 
show MAP control voltage V* (left) and the actual vibration V sr jB (right). 

In this test, MAP control criterion meets the actual vibration measurement with reasonable 
agreement, that is, the minimal vibration occurs near the secondary voltage value where 
secondary power equals the actuator aggregated loss (intersection) and the net power input 
to the plate is minimum. The phase relation of the excitations is kept very simple. Both 
secondary actuators are set to 0° and the phase of the primaries are 180°. As predicted by 
the model for low damping and single excitation, mMAP control has demonstrated small 
sensitivity to the phase of secondary actuation 


Multi Excitation with Primary Arbitrary Phases 

Unique and distinctive characteristics of mMAP control are fully revealed when arbitrary 
phasing is imposed on the primary sources. The determination of the optional control now 
relies on optimizing both the magnitude and phase. The phenomenon of phase dependency 
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is fully observed in the experimental results shown in Figures 4.20-4.21. 

The aluminum plate is excited by both primaries and only one secondary performs the control 
to illustrate the complex phasing challenge. The MAP optimal control is again indicated 
by the intersection of the secondary power to its aggregated power loss (zero net power) 
and the total input power is indicated by the sum of all actuator power. The input to the 
secondary actuator is swept while the inputs to the two primaries are kept constant. The 
vibration level is monitored by the accelerometer at the plate center. Several observations 
can be made from these test results. 

Figure 4.20(a) shows that the rnMAP computes the optimal control voltage less accurately 
with two secondary actuators than with a single secondary actuator. In Figure 4.20(a), V* 
nearly equals V S( jb while in Figure 4.19, these two voltages are slightly separated, although 
in both cases the primary phasings are simply set to zero. Figures 4.20(a)-(b) and 4.21(a) 
show the radical change of the MAP optimal control with respect to the minimal vibration 
measurement and the total power input as the primary phasing varies from 0° to 90° and the 
secondary phase is fixed to 180°. In addition, the vibration level dip near the optimal control 
is very sharp for zero primary phasing. It gradually flattens out as the primary phasing shifts 
to 90°. Similarly, the alignment of minimal input power to the minimal vibration also shifts 
but at a milder rate. 

It is interesting to note that the optimal control V* and the minimal vibration V s db flips back 
to a near perfect alignment if 45° additional phase is imposed to the secondary actuator, as 
seen in Fig. 4.21(b). Meanwhile, the vibration measurement dip also returns to this original 
sharpness as the secondary phase is optimal. Note the optimal phase 45° is half the primary 
phasing of 90°. 

The phase dependency of mMAP control is quantitatively summarized in Fig. 4.22. We use 
the ratio of the actual minimal vibration voltage V Sf jb over the MAP control voltage V,* as 
an indicator of the MAP convergence and plot it as a function of the primary phasing, with 
unity indicating perfect agreement between the two voltages. The indicator climbs to 2.5 
when the primaries are out of phase by 90°. Setting the phase to 45° brings the convergence 
indicator back to 1. These results illustrate that for mMAP control with arbitrary phasing, 
the secondary phase must be controlled and optimized while for near zero primary phasing, 
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Figure 4.20: (a) 0 degree primary phasing and (b)30 degree primary phasing. Two dashed 
vertical lines indicate the minimal vibration and MAP control, respective. V p \ = 0.5 V, 
V P 2 = 0.5 V, / = 1136 Hz, S\ phase = 180°. 

MAP control can be performed with the secondary phase neglected without a minimal control 
error. 
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Figure 4.21: (a) 90 degree primary phasing and S\ phase = 180° and (b)90 degree primary 
phasing with Si phase = 225°. V p \ = 0.5 V, V P 2 = 0.5 V, / = 1136 Hz. 

MAP Demonstration on Composite Structure 


A composite test rig was instrumented for multiple excitation control and preliminary tests 
preparation for the real time mMAP demonstration on a more realistic helicopter structure. 
The composite panel, a mock up of a helicopter fuselage skin, consists of two bays and single 
plate base panel. The composite panel consists of a panel skin with 6 plies of graphite for 
a total thickness of 0.045 inches and a Kevlar Core 0.75 inches thick and the weight of the 
panel is 14.48 lbs. The frequency response function, determined using a tap test, is used to 
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Figure 4.22: Convergence of MAP control vs primary phasing. 
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Figure 4.23: Frequency response function for the composite panel. 


extract the damping values as shown in Fig. 4.23. 

The damping value for different frequencies for the composite panel is given in Table 4.3. 
Structural damping was extracted using a multi-degree of freedom modal curve fitting tech- 
nique within the IDEAS software. 

Four piezoelectric patch actuator/sensors are mounted on symmetrically about the center of 
the panel as primary and secondary actuators as shown in Fig. 4.24. The panel is hung from 
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Table 4.3: Damping versus frequency for the composite panel 


Freq(Hz) 

7 1% 

293 

3.0 

539 

2.6 

874 

4.0 

1068 

8.9 

1473 

4.8 




Figure 4.24: Composite skin panel with two primary and two secondary piezoelectric patch 
actuator/sensors mounted. 

two bungee cords to simulate a free-free boundary conditions. Actuators are mounted on 
the panel using an improved electroding method in which the copper foil is directly soldered 
to the PZT electrode. The actuators are driven by audio amplifiers. 

Compared to the aluminum panel tested previously, two distinctive power consumption 
features are observed. First, the composite panel shows heavier structural damping as high 
as 10% so that the dynamic range of the electric power (on resonance vs off resonance) of 
each actuator around a resonance is significantly reduced (much smaller resonance power 
peaks). Second, the off- resonance power changes by a factor of four from 400 Hz to 1000 
Hz, as seen in Fig. 4.25. While both panel damping and the piezoelectric dielectric loss 
are responsible for this frequency dependence, the former is believed to be dominant as the 
piezoelectric loss varies by less than 10% (as determined experimentally) in this frequency 
range. The power spectrum also shows there is very little identifiable modal activity beyond 
1000 Hz due to heavy damping. 
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Figure 4.25: Electric power measurement of four actuator-sensor on the composite panel at 
10 V. 

For MAP control, the high structural damping constitutes a major challenge. Effective sens- 
ing and control of vibration by the secondary actuators relies on the structure’s transfer 
mobility which is a function of material damping and structural configuration. The high 
damping weakens the communication between the primary and secondary. Figure 4.26 ana- 
lytically shows the effect of the high structural damping on the MAP control. As the damping 
increases, the secondary actuator power approaches its aggregated loss in shape and value 
and the intersection of them becomes indistinct. This could result in a large uncertainty 
in experimental determination of the MAP control point when both the control circuit and 
aggregated loss have a limited signal to noise ratio. 

Nonetheless, to validate the controllability on a high damping structure, we have performed 
a preliminary MAP test on the composite panel using the single primary and single sec- 
ondary configuration. Figure 4.27 shows that the vibration reduction is comparable to the 
low damping case (-12 dB), however, the result agrees with the simulation results. The 
aggregated loss and secondary power curves run nearly parallel to each other in the figure 
and the experimental noise creates multiple intersections between these curves. Therefore, 
detecting the optimal control voltage becomes a challenge. 
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(a) (b) 

Figure 4.26: Effect of structural damping on MAP control: (a) low damping 1% (b) com- 
posite high damping 10%. 


To enhance MAP implementation, two techniques are being sought. One is to increase the 
primary excitation strength and to force more communication between the primary and 
secondary. The second technique is to aggressively filter the aggregated loss in order to 
obtain clean intersections. Figure 4.28 shows an effort to enhance the MAP determination 
with heavy damping loss. When the primary excitation is doubled, the secondary power P s2 
is now more negative and the two power curves are more separated compared to the Vp=5 
case. However, because of the instrument voltage range limitation, the new intersection 
between the secondary power and secondary loss cannot be observed. We anticipate that 
the two curves intersect near 15 V. 


130 



Figure 4.27: MAP control on high damping panel (V p = 5 V, center accel). 



Figure 4.28: Effect of augmented primary excitation on MAP control V p — 10 V, Center 
Accel. P p i is out of range. Power scale is set to compare with V p = 5 V case. 

4.4 Summary 


We summarize the key accomplishments presented in this chapter: 


• The theory of MAP was extended for multiple primary excitations controlled by mul- 
tiple secondary excitations. 
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• The effect of primary phasing on MAP control was investigated. It was shown that 
when the primary sources have zero relative phasing, the real secondary power is zero 
at optimal control. However, when the primary sources have arbitrary relative phasing, 
the secondary real power need not be zero at optimal conditions. Symmetry conditions 
under which the secondary real power is guaranteed to be zero at optimal control were 
derived. However, satisfying the symmetry constraints on the secondary locations with 
respect to the primary sources may be difficult in practice. 

• An electro-mechanical model for a simply supported plate excited by multiple primary 
and multiple secondary PZT patch actuators was developed. 

• MAP control with multiple primary and multiple secondary PZT actuators was demon- 
strated experimentally on a composite panel 
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Chapter 5 

MAP Feasibility for Rotorcraft 
Vibration Control 


In chapter §4, we developed the theory of MAP for multiple primary sources controlled by 
multiple secondary sources. It was shown that when the phasing between the primary sources 
is zero, then the secondary real power is zero at optimal control. However, when the primary 
sources have arbitrary relative phasing, then the secondary real power is not zero. This poses 
a challenge as MAP control is hard to implement if the secondary real power is not zero 
at optimal control. Also, when the primary sources have arbitrary phasing, computing the 
optimal secondary phase is not straight forward. Finally, extending MAP control to multiple 
secondary sources requires a feasible approach to find the optimal solution especially when 
the secondary sources are numerous. The main focus of this chapter is to assess the feasibility 
of MAP control for realistic structures having a distributed primary excitation with arbitrary 
relative phasing and possibly multiple frequencies controlled by a large number of secondary 
sources. 
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5.1 MAP Solution for Arbitrary Primary Phasing 


In this section, we present new theoretical results for MAP control implementation for mul- 
tiple primary sources with completely arbitrary phasing. In §5.1.1, we present the theory 
for multiple primary sources and one secondary source. We show that placing the secondary 
actuator at the intersection of the node lines of the neighboring modes, that is, the modes 
before and after the mode of interest, guarantees that the secondary real power is zero at 
optimal condition. We also show that the optimal secondary phase and amplitude can be ex- 
actly calculated without any knowledge about the primary sources and phasing at resonance 
and off-resonance frequencies. In §5.1.2, we present numerical simulations to validate the 
theoretical findings for resonance and off-resonance cases. In §5.1.3, we describe an iterative 
algorithm to find the optimal secondary sources and prove asymptotic convergence of the 
algorithm to the optimal solution. However, this approach is not feasible for large number 
of secondary sources. We present an alternative approach to find the optimal solution for 
large number of secondary sources. 


5.1.1 MAP Solution for Multiple Primary and Single Secondary 

Consider a two-dimensional structure that is excited by N primary sources (point forces) 
Pi,P2, ■ ■ ■ ,Pn ■ The primary sources are assumed to have arbitrary phasing and given by 

Fp=[F P i F p 2 ... F pN ] = [f pl e^ 1 f p2 e^ ... f pN e?**»] 

where f Pi and <p Pi are, respectively, the amplitude and phase of the i th primary source. 
Assume that the structure is controlled by one secondary source given by 

F s = f s e j *° 

where f s and (p s are the amplitude and phase of the secondary excitation. Note that the 
optimal secondary source is 

f; = /,V« = -R-}R, r F p (5.1) 
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where /* and </>* are the optimal secondary amplitude and phase, R ss is the real part of the 
point mobility and R sp is the real part of the transfer mobility M sp given by 

R sp [R SP1 R , sp2 ’ ’ ’ -h > . S pv j • (5.2) 

The MAP condition requires that the optimal secondary solution (/*, </>*) lies on the P s = 0 
curve where P s is the real secondary input power. In general, for arbitrary primary phasing, 
the optimal solution does not lie on the P s = 0 curve. However, we show that optimizing the 
location of the secondary actuator guarantees that the optimal secondary solution lies on the 
P s = 0 curve. This is a significant step towards implementation of MAP control technique 
for realistic excitations. 


Consider the transfer mobility between the secondary source and the 7 th primary source 


M spi 


OO 

^(riuu 2 n + ju{u 2 n -u 2 )) 

n = 1 


40 n (x s , y s )(j) n {xp il y Pi ) 
m [(< 4 - (u 2 ) 2 + uiy 2 } 


(5.3) 


where n represents the mode index ( 77 , 1 , 77 . 2 ), c 0 is the excitation frequency, c o n is the 77 th 
resonant frequency, r/ is the damping, cj) n is the 77 th mode shape, (x Pi ,y Pi ) is the location of 
the 7 th primary and (x s ,y s ) is the secondary location. When c 0 is close to Uk, then most of 
the modal contribution can be assumed to come from the resonant frequencies c Ok, ^k - 1 and 
oJk+i- Equation (5.3) can be written as 


fc+i 

M spi ~ _ w2 )) 

n=k — 1 


4 (f)n{%s-) Vs) ( t > n{%pii Vpi ) 
777 [(cu 2 - OJ 2 ) 2 + OJ^r] 2 } 


(5.4) 


Now, if we choose the location of the secondary source such that it lies on the intersection 
of the nodal lines of the mode shapes corresponding to 0Jk-i and 0Jk+i, then we have 


4> k -i(x s ,y s ) = 0 & (j) k+ i(x s ,y s ) = 0 


Using (5.5) in (5.4), we get 


M spi = (r]ujuj 2 k + ju(uj i 

= ( a + jb)5 m 


w 


'■)) 


Vs) ( Pk(%pi 1 Vpi) 
m [(u>l - w 2 ) 2 + UpT] 2 


(5.5) 


(5.6) 

(5.7) 
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where 


a = r]ujuj 2 k , b = u(uj 2 k - uj 2 ), 5 spi = 


4 (ftkyXsi y s )4>k{z P i, Vpi, 


m IM - ^ 2 ) 2 + utv 2 } 

Next, we show that given (5.7), the optimal secondary solution lies on the P s 
The real secondary power can be expressed as 

P s = ^Re (F^MspFp) + ^F”R ss F s 

where using (5.7), M sp can be expressed as 

■M-sp (®T jb~j $sp2 ■ ■ ■ (u T jfySgp 

which gives 

Rsp ® [ b | & sp2 ' ' ' $s P v ] 

Using (5.1) in (5.9), we get the secondary power at optimal condition as 


P* — 
s 2 R 


1 Re ( F p R ps M sp F p ) + -^—F } H R ps R sp F p 

Zibc 


Using (5.10) and (5.11), it can be seen that 

R ps M sp a ( o T jb)5 ps 5 sp 

is symmetric. Since R ps M sp is symmetric, (5.12) reduces to 


p* _ 


2 R„ 


(. F p R ps R sp F p ) 


f h r r f 


2 R„ 


= 0 


(5.8) 
0 curve. 

(5.9) 

(5.10) 

(5.11) 

(5.12) 

(5.13) 


(5.14) 

(5.15) 


Thus, the optimal secondary solution lies on the P s = 0 curve. Now, if we know the optimal 
secondary phase, then the optimal secondary amplitude can be easily calculated by sweeping 
the secondary amplitude. Next, we show how the secondary optimal phase can be computed. 


Note that M sp , as given in (5.10), can be expressed as 


M sp = [ft ft • • • ftv] 


(5.16) 
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where fa = mag((a + jb)S spi )) and 



_i b _i (cUfc - w 2 ) 

-0 = tan - = tan 1 5 — 

a r]U} z k 

(5.17) 

which gives 

R sp =[f3 1 /? 2 ••• /?iv] cos 0 

(5.18) 

Using (5.1) and (5.18), the optimal secondary source can be written as 



N 

f* e i<t>t — cos -0 E 
2=1 

(5.19) 

where 

A/p* 

P2 t~> 

^SS 

(5.20) 

Rearranging, we get 

N 

f* = cos-0 E 
2=1 

(5.21) 


Since /* is real, setting the imaginary part to zero provides a condition that the optimal 
phase has to satisfy 

N 

E d* sin (<^ - o = o ( 5 - 22 ) 

The real secondary power can be expressed as 

1 1 N 

P s = -Rssfs + 2 E fsfofpi cos (<Ai + ^ - 0 S ) (5-23) 

2=1 

Setting P s = 0 and expressing f s in terms of we get the equation for the P s = 0 curve 

fa{<j>a) = E-^^COS (<l>i + ~ </>a) 

-i -Li'SS 

2=1 

N 

= E^ c °s(^ + ^-^) (5.24) 

2=1 
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(5.25) 


Differentiating (5.24) with respect to 0 S gives 


df s 


N 


d<t>, tr 

If (j) s is a solution to the equation 


Y fJ-i sin (fa + 0 
dfs 


= 0 


then, (5.25) gives 


N 


Y ^ sin (^ + ip - 00 = o 


i=i 


Comparing (5.22) and (5.27), the optimal secondary phase is given by 

0* = 0 S ~ if or 0* = 7T + (j) s - 0 
A special case is oj = u>k (resonance). In this case, 


0 = tan 

which gives the optimal secondary phase as 


-i H - u 2 ) 




= 0 


's — Ys 


or 0 S = 7T + 


(5.26) 

(5.27) 

(5.28) 

(5.29) 

(5.30) 


The new MAP algorithm for multiple primary sources and one secondary source is summa- 
rized: 


• Suppose that u ~ uik- Place the secondary source at the nodal intersection of the mode 
shapes for uJk-i and cok+i- 

• Find the phase 0 S where the P s = 0 curve is a minimum or a maximum. 

• Calculate 0 according to (5.17). 

• Calculate the optimal secondary phase 0* according to (5.28). 

• Set 0 S = 0* and sweep f s till P s = 0 to give the optimal /*. 
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Figure 5.1: Primary source locations shown with corresponding amplitude and phase. Lo- 
cation A corresponds to the optimal secondary location while location B is not optimal 

5.1.2 Simulation Results 

We consider three primary point excitations with arbitrary phasing and one secondary point 
excitation acting on the simply-supported aluminum plate previously considered in other 
simulation studies. The primary sources with the relative phasing are shown in Fig 5.1. 
We assume that the excitation frequency is close to the resonant frequency of oj k = 826 
Hz. The corresponding mode number is (1,3). Note that the neighboring frequencies are 
ojk-i = 724 Hz with mode number (3,1) and cu k +i = 976 Hz with mode number (3,2). 
Therefore, according to the theory, the optimal secondary location is at the intersection of the 
nodal lines for modes (3,1) and (3,2) as shown by location A in Figure 1. Location B shows 
a non-optimal secondary location. We consider three cases: (1) oj = 826 Hz (resonance), (2) 
(jj = 800 Hz (off-resonance and left of 826 Hz) and co = 850 Hz (off-resonance and right of 
826 Hz). 

Figure 5.2 shows the simulation for resonance case where P s = 0 line is plotted in the f s — </> s 
space. The contour lines for the total input power are superimposed. In Figure 5.2(a), it can 
be seen that the turning point of P s = 0 curve (equivalent to minimum or maximum points 
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Figure 5.2: Resonance case u = 826 Hz. The P s = 0 line is shown in dark. The contour 
lines for the total input power are superimposed, (a) Secondary located at A (b) Secondary 
located at B 

when phase is plotted on the x-axis) corresponds to the optimal secondary shown in red, as 
suggested by the theoretical results. As expected, the total input power is minimized at the 
optimal condition. Figure 5.2(b) shows the simulation for the non-optimal location. It can 
be seen that the optimal solution (shown in red) does not lie on the P s = 0 curve. 

Figures (5.3) and (5.4) show the simulation for the off-resonance cases. The calculated 
optimal phase is indicated by the blue line. The intersection of the blue line with the P s = 0 
curve gives the estimated optimal secondary source which is close to the true optimal given 
by the red dot. 

Next, we try to understand the trade-offs between total power attenuation and MAP im- 
plementation based on optimal secondary location with respect to the mode shapes. We 
consider a simply supported plate that is excited by 8 primary sources (point forces) along 
the edges as shown in Fig 5.5. The primary locations are indicated by red dots. The relative 
strength and phasing between the primary sources is arbitrary. For the simulations, we chose 
the primary sources to have the same amplitude but different phasing. Specifically, the phas- 
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Figure 5.3: Off-resonance case u = 850 Hz. The P s = 0 line is shown in dark. The contour 
lines for the total input power are superimposed. 

ing was p = [0, 45, 30, 20, -30, -30, -60, 120] degrees with the ordering starting from bottom 
left and moving counter-clockwise. We consider off-resonance excitation frequency of 800 Hz 
which is close to the resonant frequency of 826 Hz. The secondary source location is varied in 
the interior of the plate and the optimal control is applied at each location assuming that the 
primary sources are known. The total power attenuation P re d = ( P nc — P c ) / P nc , where P nc 
is the total input power without any control and P c is the input power with optimal control, 
is plotted. Note that P red =l for no control and P re d= 0 for complete control. We compare 
the power attenuation at two different locations. Location A is optimal with respect to the 
mode shapes adjacent to 826 Hz. Location B is the plate center and hot spot for 826 Hz. It 
can be seen that the power attenuation at A is close to the maximum attenuation you can 
get. If one uses location B for the secondary, then we will get good power attenuation in the 
ideal case where the primary sources are known. However, for MAP control we rely only on 
the secondary power for control. Fig 5.6 shows the variation of the real secondary power at 
optimal control. 
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Figure 5.4: Off-resonance case to = 800 Hz. The P s = 0 line is shown in dark. The contour 
lines for the total input power are superimposed. 


It can be seen that the secondary real power at A is zero, as predicted by the theory, and 
hence MAP control can be implemented. However, at location B, the secondary real power 
is not close to zero and therefore MAP control is not suitable even though the attenuation 
is very good at B. 


5.1.3 MAP Solution for Multiple Primary and Multiple Secondary 
Sources 

In §5.1.1, we presented an algorithm to find the optimal secondary source when multiple 
primary sources with arbitrary phasing are present. Given that we know the optimal solution 
for one secondary source, we present two approaches to find the MAP solution for multiple 
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Figure 5.5: Simply supported plate excited by 8 primary sources with excitation frequency 
equal to 800 Hz. The total power attenuation is shown as a function of the secondary location. 
Location A shows the optimal location of secondary with respect to the neighboring mode 
shapes. Location B is plate center and hot spot for the mode at 826 Hz. 

secondary sources. We first present an iterative algorithm to find the optimal secondary 
sources when multiple secondary sources are used for MAP control. We show that the 
iterative algorithm is asymptotically stable and converges to the optimal value. The iterative 
solution, however, is not feasible for large number of secondary sources. We present an 
alternative approach that is feasible for large number of secondary sources and more useful 
in practice. 
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Figure 5.6: Simply supported plate excited by 8 primary sources with excitation frequency 
equal to 800 Hz. The real secondary power is shown as a function of the secondary location. 

Iterative Algorithm 

Without loss of generality, consider N primary sources (point forces) Pi,P2, ■ ■ ■ ,Pn with 
arbitrary phasing 

Fp = [F p i F p2 ... F pN \ = \f pl e^ 1 f p2 e- *** ... f pN e^} (5.31) 

For purposes of illustration, we consider two secondary sources si and s 2 given by 

F?=[F s1 F s2 \ = [f sl e^ sl f s2 e^} (5.32) 

Recall that the optimal secondary sources are given by 

F: = -R-'RspFp (5.33) 


144 


where the point and cross-mobility R ss and R sp are given by 



(5.34) 


and 




(5.35) 


The iterative algorithm is as follows: 


1. Assuming only si to be active, find the optimal F° . 

2. Fix Si to be F ' s ° and hnd the optimal s 2 , -F s ° 2 , by treating si as a primary source. 

3. Fix S 2 to be F s ° 2 and hnd the optimal Si, F ^ , by treating S 2 as a primary source. 

4. Fix si to be F^ and hnd the optimal S 2 , i 7 ^, by treating si as a primary source. 

5. Repeat the iteration procedure till \F^ +1 — F*| < e where i — 1,2 and e is some 

convergence metric. 

The iterative procedure can be represented as a map: 


for i — 1,2. Next, we show that the iterative map given in Eq. 5.36 is asymptotically stable 
and converges to the true optimal value given in Eq 5.33. 


Stability of the Iterative Algorithm 

Assuming si only, the optimal value of Si is given by 


pk pk + 1 


(5.36) 




(5.37) 
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(5.38) 


Treating si as a primary source equal to F s ° , we get the optimal s 2 as 


J? = -§^F„ - = o 2 + &F° 


s 2 


JD ~ V JD ^ Si 
ri S2S2 ri S2S2 


Si 


Treating s 2 as a primary source equal to F s ° 2 , we get the optimal si as 


U = a '~R 


SlS 2 j-iO 


S 2 


SlSl 

= «i + 


S 2 


Using (5.38) in (5.39), we get 


Fg — + /3 i « 2 + PifaF 


Similarly, Fi can be expressed as 


Fg 2 — a 2 + /3 2 q?i + /?i/ 3 2 F s ° 2 
From (5.40) and (5.41), in general, we have the following iterative map: 

Fg +1 = + /3iCt: 2 + PifoFg 


F l +1 = a 2 + foaii + /3 1( 5 2 F 


(5.39) 


(5.40) 


(5.41) 


(5.42) 


Next, we show that the steady state solution of the iterative map given in (5.42) corresponds 
to the true optimal solution. From (5.42), the steady state solution satisfies 


which gives 


F ’ s * — oil + P 1012 + @102 F* ± 

F* 2 = ct 2 + /l 2 «i + fii/3 2 F* 2 


= rv^s( a2+ftai) 


(5.43) 


(5.44) 
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Using expressions for cki,0! 2 ,/3i and /3 2 , we get 


\K] 


i 

Rs 2 S 2 

Rs\S 2 



Rs 1 siRs 2 s 2 — R‘sis 2 

R-sis 2 

Rs 1S1 


= -R-jRs P F r (5.45) 


which shows that the equilibrium point of the iterative map is the optimal secondary source 
distribution. Next, we show that the iterative scheme is asymptotically stable. Consider the 
map for F Sl 

F s fc + 1 = a 1 + /3 1 a 2 + /Wj 

The map is linear and the eigenvalue is A = [5 i/ 3 2 . Asymptotic stability is guaranteed if 
A < 1. Therefore, we need 

aft = < i 

-Ll'SiSi - r *'S2S2 

Since R ss is positive definite, we have 

A, S i. sl R,s 2 S2 - Rsi S 2 > 0 =► p 'p 2 < 1 

-“-Si Si t*'S2S2 

which guarantees asymptotic stability. 


Alternative Solution for Multiple Secondary Sources 


The iterative approach presented in the previous section is not feasible for large number of 
secondary sources. For large number of secondary sources, we need a more practical solution. 
In what follows, we present an alternative approach that is feasible when the secondary 
sources are large in number. Since the location of the secondary sources in known, we 
assume we can compute R ss experimentally before running MAP control. Recall that the 
optimal secondary sources are 

F: = —RJs RspF p (5.46) 

Note that we assume R ss is known. However, R sp F p is not known. First, compute the 
optimal secondary sources one at a time. That is, 
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1 


(5.47) 


F Si = - n Rs :P F p V j = 1 , 2, . . . , m 


R* 


Writing (5.47) compactly, we get 


F s diag(J? ss ) R sp Fp 


where 

Fj = [F Sl F S2 ... F Sm } 

Equation (5.48) can be re-written as 

RspFp diag( R ss ^jF s 

Using (5.49) in (5.46), we get the optimal secondary sources as 

F* = R~gdiag(R ss )F s 


(5.48) 


(5.49) 


(5.50) 


Thus, for m secondary sources, one needs to run MAP test m times to calculate the optimal 
secondary sources. 


5.1.4 FEM Model 


A FEM model was built to verify MAP control authority and global vibration reduction 
and also guide the experimental test rig design and validation. The design consists of a 
thin aluminum plate riveted to rigid aluminum angles and is representative of a realistic 
airframe structure. Point force or line moment excitations were applied to the angles as 
primary excitation. The force amplitude and phase of a secondary point force excitation 
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Secondary 

source 


Figure 5.7: Finite element model of test rig. 


were swept until the non trivial MAP point was found. The location of the secondary source 
was not optimized. A picture of the finite element model is presented in Fig. 5.7. Numerical 
boundary conditions were free-free. 

Flexural kinetic energies for a single primary point force excitation for both the system and 
the plate subcomponent are presented in Fig. 5.8. Frequency responses exhibit many modes 
in the 500 to 2500 Hz interval. The system’s flexural kinetic energy with and without MAP 
control around three of the structure’s resonant frequencies (937, 1012, and 1215 Hz) are 
presented in Fig. 5. 9. MAP produced over 8 dB of attenuation at 937 and 1215 Hz, and 
no control at 1012 Hz. This lack of control for that particular resonance was due to the 
secondary excitation being located on a nodal line of the underlying structural mode. 

The effect of the MAP control on the radiated sound held is presented in Fig. 5.10. When 
modes could be controlled, the radiated sound power was reduced by 9.5 and 8.5 dB at 937 
Hz and 1215 Hz, respectively. Figures 5.8 and 5.9 show that there is a direct relationship 
between structural vibration and interior noise. 
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Figure 5.8: Rig flexural kinetic energy 
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Figure 5.9: Flexural Kinetic Energy with and without MAP control. 

5.2 Power Model for a Generic Actuator 


The power model of a generic PZT actuator was developed in §3.4. Specifically, we relate 
the net electrical power in the actuator, the quantity that is measured, to the net mechanical 
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Figure 5.10: Radiated sound power with and without MAP control. 


power, the quantity that is controlled, through the dielectric loss using the following equation: 

1 oj a 2 


Pe 2 h 


|V| 2 e 33 6 + P 


(5.51) 


Equation (5.51) shows that the net input electrical power P e is equal to the sum of the 
dielectric loss and the net input mechanical power P. Note that even though Eq. (5.51) 
was derived for a PZT actuator, a similar relation can be derived for other actuators, for 
example, proof-mass actuators. 


5.3 Test Rig and Plate Characterization 


The main experimental efforts this year include: 1) verification of global power minimization 
using MAP, 2) implementation of MAP control on a realistic and sophisticated aircraft 
structure, and 3) experimental validation of the new MAP theory developed in §5.1. 
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Figure 5.11: Test rig with a shaker and a PZT sensor /actuator pair. 


A test apparatus was constructed to mimic a representative rotorcraft skin structure with 
more realistic vibration load paths from bulkhead to skin. An aluminum plate, 12 x 12 x 1/8 
inches, was riveted by 28 aircraft grade rivets to a welded lxl inch aluminum angle frame. 
This design introduces contact surface and Coulomb damping to the test rig. The plate 
assembly is suspended by 4 bungee cords to mimic a free-free boundary condition, as shown 
in Figure 5.11. Two PZT patch sensor /actuators (2.5 x 2.5 x 0.005 inches) are bonded to the 
upper left and low right corners, respectively, one acting as a primary actuator (disturbance) 
and the other acting as a secondary actuator (control). A mini-shaker is attached to the 
plate via a piano wire stinger at the aluminum side wall as an alternative primary actuator 
to mimic distributed vibration load along the edges of the side walls to the aluminum plate. 

The structural dynamic characterization of the newly constructed plate was conducted with a 
Polytec PSV-400/OFV5000 Laser Scanning Doppler Vibrometer (LSDV) as shown in Figure 
5.12. The LSDV scans the plate surface in a grid of 320 x 320 nodes, and a laser pointer dwells 
on each node for a certain period to measure average vibration velocities orthogonal to the 
plane. The plate is excited by either the shaker or the primary PZT with a broad frequency 
band noise signal. The frequency response function (FRF) of the vibration velocities with 
respect to the excitation at each node is then calculated, and mode shapes of any selected 
frequencies are extracted from the FRFs. 
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Figure 5.12: Polytec Laser Scanning Doppler Vibrometer 

FRF of Framed Plate 



Figure 5.13: FRF of the plate measured by LSVD with shaker excitation. Red circle indicates 
the frequency selected for MAP testing 


Figure 5.13 shows the measured mobility at the center of the plate for a shaker excitation at 
the edge of the side. There are over 20 strong vibration modes in the 0-3000 Hz frequency 
range. The modal cluster around 1000 Hz is of particular interest and is selected for MAP 
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control and analysis. The velocity mode shapes of 3 frequencies, 624 Hz, 999 Hz and 1052 


Figure 5.14: Velocity modal shapes for (a) 624 Hz (b) 999 Hz and (c) 1052 Hz 


5.4 Experimental Validation of MAP Control 

The validation of the global energy minimum at optimal control was conducted using a pair 
of PZT patches rather than using the shaker (as primary) and a PZT patch (secondary) due 
to some uncertainty and slight asymmetry of the mode shapes imposed by the attachment 
of the shaker. A frequency of 999 Hz (2,2 mode) was selected for MAP control. As can be 
seen from Figure 5.13, this frequency is closely coupled to other two modes. Minimization 
of global kinetic energy of the plate at optimal MAP control is verified by: 

• Measuring velocity mode shapes of the plate at various excitation levels 

• Calculating global kinetic energy of the plate 


Hz were extracted from the FRFs at all the nodes, as shown in Figure 5.14. The (2,2) mode 
of the plate at 999 Hz is selected for MAP control due to symmetry of the mode shape with 
respect to PZT actuator locations. 



(a) 


(b) 


(c) 



n 


(5.52) 


i=l j = 1 
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Figure 5.15: MAP control with a pair of PZT patch sensor/actuators. Red circle indi- 
cates MAP predicted optimal control at 14=5.94 V. Purple circle indicates optimal control 
at 14=5. 70V as predicted by the accelerometer. 

where Vij and rn i3 are the vibration velocity and the mass of plate element at location 
i and j, respectively and l x n is the LSDV scan grid size. 

• Identifying the optimal MAP control voltage V* by sweeping 14 applied to the sec- 
ondary actuator. 

• Checking the correlation between V*, minimum and minimum acceleration mea- 
surement. 


The harmonic primary excitation is kept constant at V p = 5 V with primary phase equal to 0°. 
The optimal secondary phase was estimated to be 183.8° with respect to the primary. This 
is close to the true optimal of 180° (or 0°). The secondary voltage 14 is swept from 0 to 10 
volts while its phase is kept fixed at 183.8°. The secondary power is acquired and compared 
to the predetermined aggregated power loss. Meanwhile, the acceleration of the plate center 
is monitored by the accelerometer. Figure 5.15 shows that the MAP criterion (intersection 
of the baseline with the secondary real electric power) predicts the optimal control to be 
14 = 5.94 V while the acceleration reading indicates maximum vibration reduction(from -26 
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Figure 5.16: Normalized total kinetic power of the plate around optimal MAP control. Blue 
circle indicates global kinetic power of the plate at V s = 5.60 V, and the normalized total 
energy is 28.1 at no control ( V s = 0). 

dB to -42 dB) at V s = 5.70 V. Figure 5.16, from the LSDV measurements, shows the total 
kinetic energy of the plate calculated from the mode shapes at different secondary voltages 
and normalized with respect to the total power at optimal control, that is, at V s = 5.94 V 
and shows that the global power minimum occurs when V s = 5.60 V. It can be seen that 
the acceleration measurement predicts the global energy minimum very well even though it 
is slightly off the hot spot of the mode shape. 


5.4.1 Change in Mode Shapes with MAP Optimal Control 


The effect of MAP control on mode shapes of the plate at 999 Hz is shown in Figure 5.17. The 
MAP control significantly reduces the vibration magnitude in the entire plate. It was also 
observed that the PZT actuators placed on the hot spots of the mode shape do not change 
the mode shape of a resonance frequency while significantly reducing the amplitude of the 
vibration. This test has demonstrated that optimal MAP control corresponds to minimum 
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(a) No MAP Control (b) With MAP Control 


Figure 5.17: Velocity modal shapes in mm/s for no-control condition and with optimal 
MAP control (V s = 5.5V). The low left square is the primary (excitation) and upper right 
the secondary (control). 

global kinematic energy of the plate. The test also demonstrates that MAP control is capable 
of reducing the vibration of a plate with complexity in structure and damping mechanism 
effectively. 


5.4.2 Validation of New MAP Algorithm 

The main idea behind the new algorithm is to optimally locate the secondary actuator 
such that the cross mobility terms can be approximated as a single mode. This guarantees 
that MAP solution holds, that is, the optimal solution lies on the P s = 0 curve for any 
arbitrary phasing of the primary actuators. The optimal secondary amplitude and phase 
can be determined for resonance and off-resonance cases. We used the shaker as the primary 
excitation and one of the PZT patches as the secondary actuator as shown in Fig. 5.11, §5.3. 
The excitation frequency was kept at 994 Hz. Figure 5.18 shows the experimental mode 
shape corresponding to 994 Hz. The secondary PZT actuator sq is shown in bottom left and 
is located close to the hot spot and is optimal for 994 Hz. 

The first step in the new MAP algorithm is to find the optimal secondary phase. For resonant 
excitation, the optimal phase is given by the secondary phase at which the secondary power 


157 



Figure 5.18: Experimental mode shape corresponding to resonant frequency 994 Hz 



Figure 5.19: Plot showing implementation of optimal phase search. The optimal secondary 
phase is the phase at which the secondary phase is minimum or maximum. 


is either a minimum or maximum for a given voltage. The primary excitation was 5 V and 
the phase was set to zero. With the secondary voltage equal to 3 V, the secondary phase was 
varied to find the phase for which the secondary power is a minimum or maximum. Figure 


Vibration (dB) 



5.19 shows the variation of the secondary power P s with respect to the secondary phase. It 
can be seen that the estimated optimal phase is —31° where P s is minimum, indicated by 
the red circle in Figure 5.19. Figure 5.19 also shows that the accelerometer reading, located 
close to the plate center, indicates that the vibration minimum occurs at —26° (purple circle) 
which is close to the estimated value of —31°. 



Figure 5.20: MAP control with secondary phase = —31° 


The next step is to fix the secondary phase at the optimal phase and sweep the secondary 
voltage to get the optimal voltage, indicated by intersection of the secondary power and the 
baseline loss curve. Figure 5.20 shows MAP performance when the secondary phase was set 
to —31°. It can be seen from the accelerometer reading that a vibration reduction of 20 dB 
is achieved. When the secondary phase was set to —26°, the MAP performance was greatly 
enhanced leading to a vibration reduction of 44 dB, see Figure 5.21. 

In summary, the new MAP algorithm provides a novel way to find the optimal secondary 
phase and amplitude without any knowledge about the primary sources. 
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Figure 5.21: MAP control with secondary phase = —26° 

5.5 Rotorcraft Airframe Mockup 


Modal analysis of the riveted structure revealed asymmetrical mode shapes and split reso- 
nance peaks that likely were resulting from asymmetrical boundary conditions due to the 
rivets and small variations in this fabricated structure. These experimental modes shapes 
could not be easily predicted using finite element models. To position secondary actuators 
for MAP, it is best to know the structural mode shapes to guarantee optimal control. Hence, 
to minimize the uncertainty in the mode shapes, we constructed a test structure without 
rivets but of similar geometry than the original panel. This new structure could be easily 
modeled and had predictable mode shapes and good separation of resonance frequencies. 
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Figure 5.22: Test Plate Dimensions 


5.5.1 Test Rig 

The new test plate was fabricated through a milling process to enable seamless integration 
of angle frames to the plate. Detailed dimensions of the test plate are shown in Fig. 5.22. 
The assembly is suspended in air by 4 bungee cords to mimic free-free boundary conditions, 
as shown in Fig. 5.23. A mini-shaker is attached to the plate via a piano wire stinger at the 
aluminum side wall to mimic distributed vibration load along the edge of the plate. 
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(a) 


(b) 


Figure 5.23: Picture of Test Rig showing (a) Front Face and (b) Back Face 


5.5.2 Comparison between experimental mode shapes and FEM 
Analysis 


The velocity frequency response function (FRF) of the plate is acquired by Laser Scanning 
Doppler Vibrometer (LSDV) and the plate is excited with broad band noise by the mini- 
shaker. Figure 5.24(a) shows FRF of the plate measured by LSDV and shaker excitation. A 
finite element analysis (FEA) of the plate was done to compare the frequency response and 
mode shapes. The frequency response from FEA is shown in Fig. 5.24(b). The experimental 
modal frequencies were within 1% of FEA results for frequencies within 500-2000 Hz range. 
We also compared the mode shapes for different resonant frequencies. For example, Fig. 
5.25(a) shows the experimental mode shape, a (2,2) mode, for 947 Hz. Figure 5.25(b) shows 
that mode shape obtained from the FEA model at 940 Hz which is in excellent agreement 
with the experimental mode shape. Such good agreement is also evident at other resonant 
frequencies. For example, figures 5.26 and 5.27 compare the experimental mode shape with 
the mode shape calculated using FEA for 850 Hz and 1073 Hz respectively. 
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Figure 5.24: (a) FRF of the plate measured by LSVD and shaker excitation (b) Total plate 
KE per N of force, from FEA 



(a) (b) 

Figure 5.25: Mode shape comparison (a) Velocity modal shape at 947 Hz measured by LSVD 
and shaker excitation (b) Calculated vibration mode shape at 940 Hz 

5.6 Effectiveness of MAP Control 


In this section, we present experimental results that demonstrate the effectiveness of MAP 
control for single and multi-frequency excitation. We use the shaker as the primary excitation 
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(a) 


(b) 


Figure 5.26: Mode shape comparison (a) Velocity modal shape at 850 Hz measured by LSVD 
and shaker excitation (b) Calculated vibration mode shape at 850 Hz 


fK n \ 


(a) 



Figure 5.27: Mode shape comparison (a) Velocity modal shape at 1073 Hz measured by 
LSVD and shaker excitation (b) Calculated vibration mode shape at 1073 Hz 


source and a PZT patch, located at the plate center, as the secondary excitation. The 
primary excitation has two resonant frequencies, 850 Hz and 1073 Hz. Figure 5.26 shows 
that the PZT location is optimal with respect with respect to the mode at 850 Hz whereas 
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Table 5.1: Primary Excitation Parameters 


^pll 

850 Hz 

^pl2 

1073 Hz 


2.5 V 

V P 1 2 

4 V 

4>pii 

30° 

f- 

to 

0° 



Figure 5.28: Voltage sweep showing the optimal secondary voltage at 850 Hz while control 
for 1073 Hz is off. 

figure 5.27 shows that the PZT location, though good, is not the best location for the mode 
at 1073 Hz. As we discuss later, this is a good test of the effectiveness of MAP for sub- 
optimal locations of the PZT actuator. The primary excitation parameters are shown in 
Table 5.1. An accelerometer located at the plate center measured the local vibration levels 
whereas the global vibration was measured using the laser vibrometer. 


850 Hz control with 1073 Hz not controlled 

We use a frequency-by- frequency approach for controlling multi-frequency excitation. First, 
the optimal secondary voltage and phase are determined for the mode at 850 Hz while the 
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mode at 1073 Hz is not controlled. The optimal secondary phase corresponding to 850 Hz was 
found to (j) s ii = 25.09°. With the secondary phase equal to the optimal phase, the voltage is 
varied till the secondary real power intersects the baseline curve to give the optimal voltage. 
See Fig. 5.28. It can be seen that at the optimal voltage, the local vibration is a minimum 
resulting in a 30 dB reduction. 


1073 Hz control with 850 Hz controlled 

With the control for 850 Hz on, the mode at 1073 Hz is controlled next. In this case, the 
optimal phase was determined to be 0 s i 2 = —6.198°. With the secondary phase equal to the 
optimal phase, the voltage was varied to find the optimal secondary voltage at 1073 Hz as 
seen in Fig. 5.29. The local vibration level is minimum at optimal secondary voltage and 
the local vibration reduction is close to 20 dB. 



Figure 5.29: Voltage sweep showing the optimal secondary voltage at 1073 Hz while control 
for 850 Hz is active. 
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850 Hz control with 1073 Hz controlled 


Finally, we control the 850 Hz excitation again while 1,073 Hz is under control. The optimal 
phase was estimated to be (j ) s n = 26.010° with local vibration reduction of 38.31 dB, which 
is 8.31 dB more than the single frequency case. Note that the local vibration reduction 
improved as a result of better optimal phase estimation. See Fig. 5.30. The vibration 
reduction for 1073 Hz was unaffected. 



Figure 5.30: Voltage sweep showing the optimal secondary voltage at 850 Hz while control 
for 1073 Hz is on. 


Global Vibration Reduction 

Figures 5.29-5.30 show that MAP control results in significant local vibration reduction 
at both 850 Hz and 1073 Hz. However, it is important to verify whether we get similar 
performance globally. The global vibration levels were measured using the laser vibrometer. 
Figure 5.31 show the global vibration in dB with and without control. Figure 5.32 show 
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Figure 5.31: Plot showing global vibration levels with and without MAP control 


the global vibration reduction. It can be seen that vibration reduction is around 22-23 dB 
for 850 Hz excitation whereas the 1073 Hz excitation was reduced by 11-12 dB. Though the 
global vibration reduction is not as high as the local reduction, as expected, it is still very 
significant. 


Global Energy Reduction 



850 Hz excitation 1073 Hz excitation 850 Hz & 1073 Hz 

excitations 


Figure 5.32: Global vibration reduction with MAP control. 
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Figure 5.33: Illustration of PZT patch locations. 


The above results demonstrate the effectiveness of MAP control for single and multi-frequency 
excitation. Also, note that dB reduction for the mode at 1073 Hz is less than the reduction at 
850 Hz. This is due to the fact that the PZT location for 1073 Hz is sub-optimal. However, 
it is important to note that a global reduction of 11-12 dB is significant and demonstrates 
effectiveness of MAP control even when the PZT location is not optimal. 


5.6.1 MAP Control with Multiple Secondary Controllers and Mul- 
tiple Frequency 

The Lab VIEW code developed for testing was modified to excite and control three frequency 
signals with arbitrary phases. The shaker was used as the primary excitation source, and 
three PZT patches shown in Fig. 5.33 were used as secondary excitation sources. Si is placed 
at the center of the plate where optimal vibration reduction for 851.1 Hz (3 by 1 mode) is 
expected. Similarly, S 2 is located at the modal peak of 1,320 Hz. P 1 and S' 3 are positioned at 
the intersection of 866, 940, and 1310 Hz based on FE analysis. The primary excitation has 
three resonant frequencies, 851.1, 945.2, and 1319.8 Hz. The primary excitation parameters 
are shown in Table 5.2. Table 5.3 shows estimated MAP parameters after phase and voltage 
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search for the MAP control. 


Table 5.2: Primary Excitation Parameters 


Fshi 

851.1 Hz 

Fsh .2 

945.2 Hz 

Fsh:i 

1319.8 Hz 

Vshi 

2.5 V 

Vsh 2 

4.0 V 

Vsh3 

3.0 V 

4>Shl 

30° 

4>Sh2 

0 ° 

4>Sh3 

60° 


Table 5.3: Estimated Control Parameters (Si, S 2 , and S 3 ) 


Fsi 

851.1 Hz 

Fs3 

945.2 Hz 

Fs2 

1319.8 Hz 

Vsi 

-0.525 V 

V S 3 

-0.59 V 

V S 2 

-1.78 V 

4>si 

26.84° 

4>S3 

135.83° 

4>S2 

58.86° 


With the control strategy established, the plate was excited, and the vibration response 
measured using the vibrometer at 99 (9 x 11 grid) locations on the test plate, phase-referenced 
to an accelerometer located at the back side of the S 3 center. The resulting mode shapes 
are captured in Fig. 5.34. The columns (left to right) summarize the different phases of 
multi-frequency MAP control: no control and three PZT control. The spatially averaged 
vibration response is depicted in Fig. 5.35. The control case reduces the global plate response 
at all three frequencies of 851.1, 945.2 and 1319.8 Hz with significant reductions. For the 
baseline case, the average vibratory responses were around 60 to 65 dB (0 dB at 1 m/s) for 
the no control case as depicted in Fig. 5.35. When control was applied to the three modes, 
the vibratory response exhibited significant reductions for each mode (19, 4.4, and 9.9 dB 
reduction, respectively). 


MAP Control with Strategically Positioned Single Secondary 


Experimental results presented here demonstrate the vibration reduction capability of a sin- 
gle secondary for multiple-frequency disturbance. The input to the shaker is kept unchanged 
(see Table 5.2), while only S 3 is used to control three frequencies. As previously described, 
S 3 is placed at the intersection of three mode shapes (851.1, 945.2 and 1319.8 Hz). Phase 
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(a) 851 Hz (Baseline) (b) 851 Hz (Control) 



(c) 945 Hz (Baseline) (d) 945 Hz (Control) 



(e) 1320 Hz (Baseline) (f) 1320 Hz (Control) 

Figure 5.34: Experimental Deflection Shapes for Three PZT Control 

and voltage search for control parameter estimation was done before global vibration mea- 
surement as shown in Table 5.4. 
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Figure 5.35: Average plate response with and without MAP control 


Table 5.4: Estimated Control Parameters (S 3 only) 


Fs31 

851.1 Hz 

FS32 

945.2 Hz 

As'3.3 

1319.8 Hz 

V S 3l 

0.755 V 

Ks32 

-0.59 V 

VS33 

0.975 V 

4>S31 

29.84° 

(pS32 

135.83° 

4>S33 

43.67° 


Figure 5.36 shows the mode shapes before and after MAP control is applied for the respective 
frequencies. The spatially averaged vibration response is depicted in Fig. 5.37. The control 
cases reduce the global plate response at all three frequencies of 851.2, 943.7 and 1204.2 Hz 
with significant reductions as summarized in Fig. 5.38. It is noticeable that even a single PZT 
(S 3 ) alone can control multiple frequencies equivalently to three PZTs to which individual 
frequencies are assigned. In such a condition, a single control placed at the intersection of 
multiple (i.e., three and more) frequency modes can be used effectively to reduce overall 
vibration levels. 


MAP for Off-Resonance Control 


The excitation frequency was set to be harmonics of 300 Hz (600, 900, and 1,200 Hz), 
simulating gearbox vibrations not associated with natural frequencies of the test rig. Only 
S 3 was used for control. The primary excitation parameters are shown in Table 5.5 and 
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(a) 851 Hz (Baseline) (b) 851 Hz (Control) 



(c) 945 Hz (Baseline) (d) 945 Hz (Control) 



(e) 1320 Hz (Baseline) (f) 1320 Hz (Control) 


Figure 5.36: Experimental Deflection Shapes for Multi-frequency control using a Single PZT 
Control 
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Figure 5.37: Average plate response with and without MAP control 



Figure 5.38: Plot showing global vibration with and without MAP control 


Table 5.6 shows estimated MAP parameters after phase and voltage search for the MAP 
control. For the off resonance control, the challenge is that the amount of possible vibration 
reduction is dependent on the controller characteristic. Figure 5.39 is the baseline curve of S 3 
PZT patch showing power directly measured from the PZT with respect to frequency. The 
red solid line is the mechanical response of the test rig which can be suppressed to baseline 
level (dotted line). At off-resonance frequencies, the amount of structural vibration is very 
close to the baseline power level, which complicates the MAP control algorithm search for 
optimal control parameters. 

Figure 5.40 shows vibration plots before and after the MAP control is turned on. A 7.5 dB 
reduction was achieved at 1200 Hz, which is close to a resonance frequency of 1204 Hz. In 
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Table 5.5: Primary Excitation Parameters 


Fshi 

600 Hz 

Fsh2 

900 Hz 

Fsh3 

1200 Hz 

Vshi 

1.0 V 

Vsh2 

6.0 V 

Vsh3 

3.0 V 

4>Shl 

30° 

4>Sh2 

0 ° 

4>Sh3 

60° 


Table 5.6: Estimated Control Parameters (S 3 only) 


Fs 31 

600 Hz 

Fs32 

900 Hz 

F.S33 

1200 Hz 

Vs31 

-5.3 V 

VS32 

0.37 V 

VS33 

0.8 V 

0S31 

O 

O 

CO 

1 

4>S32 

O 

O 

OO 

1 

4>S33 

45° 



Figure 5.39: Plot showing S 3 baseline when 6.75 V is applied. 


600 Hz and 900 Hz cases, while local mode shape changes were evident, no significant global 
reduction was made. For application of the MAP control algorithm to real structures, it 
is suggested to closely analyze the external vibration sources, and structural responses of 
the target structure. Then controller position can be determined based on a trade study to 
reduce global vibration response due to resonance and off-resonance frequencies. 
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(a) 600 Hz (Baseline) (b) 600 Hz (Control) 



(c) 900 Hz (Baseline) (d) 900 Hz (Control) 



(e) 1200 Hz (Baseline) (f) 1200 Hz (Control) 


Figure 5.40: Experimental Deflection Shapes for Off-Resonance Control using a Single PZT 
Control 
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5.7 Summary 


The key accomplishments of the MAP feasibility analysis and rotorcraft vibration reduction 
are: 


• The feasibility of the MAP control approach was assessed for multiple primary sources 
with arbitrary phasing. It was shown that optimizing the location of the secondary 
actuator with respect to the mode shapes guarantees that the secondary real power is 
zero at optimal control for resonance and off-resonance frequencies. An algorithm to 
find the optimal phase and magnitude of the secondary source was developed. 

• A stable iterative scheme was developed to find the MAP solution for multiple sec- 
ondary sources. The iterative algorithm is attractive for small number of secondary 
sources. An alternative approach was developed that is feasible for large number of 
secondary sources. 

• Designed, modeled and built a test rig to mimic rotorcraft construction. Demonstrated 
MAP control with primary excitation using shaker with multiple frequencies and mul- 
tiple secondary actuators for control. 

• This research has led to an improvement in the technology readiness level (TRL) from 
1 to 3-4. 
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Chapter 6 

Effect of MAP Vibration Control on 
Radiated Noise For a Simple Interior 
Structure 


We present some preliminary analysis showing the potential benefit of MAP control to re- 
duce noise in the passenger cabin due to the main rotor transmission excitation, the NASA 
program goal. Helicopter gearbox excitations are transmitted to the fuselage via rigid con- 
nections at the transmission well. Subsequently, the frame and skins of the fuselage transmit 
the vibrations produced by the main gearbox and radiate sound into the cabin interior. 
While representative of helicopter interior skin panel, the plate structure used for the de- 
velopment of MAP control was relatively small and had few structural modes in the 500 to 
2500 Hz range. A fully coupled vibro-acoustic finite element model of a larger, more realistic 
and modally denser structure was used to evaluate MAP potential. Although recognized as 
a rather simplistic representation of a real helicopter interior, this model is a step closer to 
evaluating MAP performance on radiated interior noise levels. Up to this point, the MAP 
development program focused on demonstrating vibration reductions that are expected to 
lead to noise reduction. 
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6.1 Model description 


An acoustic cavity of dimension 1 x 1.1 x 1.03m 3 is enclosed in an aluminum structure 
composed of 2-nnn thick panels with 20 x 20mm 2 beams in all corners. The panels are 
modeled as shells, whereas the reinforcement beams are solid elements. A schematic of the 
structure is presented in Fig. 6.1. The base of the structure was simply supported. 



Figure 6.1: Schematic of the box structure. 

A primary point force excitation P\ is applied to one panel of the structure. MAP control 
was applied using a single secondary point force. The secondary was either located on the 
same panel as the primary (51), or on an adjacent panel (S 2 ). The locations of (Si) or (S 2 ) 
were not optimized for a broad range control or single mode control, but assigned by visual 
inspection on hotspots of the forced response at 505 Hz (Pi only). 
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6.2 Results 


The performance of the controller was evaluated by using the structural kinetic energy in 
all the panels, as well as the averaged sound pressure levels (SPL) in the acoustic cavity. 
The frequency responses without control are presented in Fig. 6.2. The many peaks in the 
structural responses suggest that the modal density is already relatively high, whereas there 
is still a lot of separation between peaks in the acoustic response. 



Figure 6.2: System frequency responses for the box structure with no control. 

For this study, the structure was excited at 505 Hz, a frequency close to a resonance of the 
acoustic cavity, and also in a range where the structure is modally dense. The structural 
and acoustical fields at 505 Hz are presented in Fig. 6.3. Results when MAP control is 
applied at S i are shown in Fig. 6.4. The MAP condition, where the structural kinetic 
energy is minimized, is noted with the vertical dashed line. The total kinetic energy level 
was reduced by 6.1 dB, and the averaged SPL dropped by 13.8 dB. While the kinetic energy 
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was minimized at MAP, a greater SPL reduction (15.5 dB) could be achieved for slightly 
larger force amplitude since the acoustic resonance is not exactly at 505 Hz and MAP controls 
the underlying dominant structural mode. 



Figure 6.3: Response fields at 505 Hz (no control), (a) Structural displacement (b) Acoustic 
pressure 



Figure 6.4: The total kinetic energy and SPL for Si MAP Control 

With MAP control applied at S 2 , the total kinetic energy level was reduced by 7.0 dB, and 
the averaged SPL dropped by 13.6 dB. Thus, MAP control can be very effective at reducing 
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the injected power even when the secondary actuation is located away from the primary 
sources and in a system where energy has many paths to propagate. 

Results of the MAP control as applied on the complex structure demonstrate that MAP 
has good potential at reducing sound inside a cavity, even if the secondary forces are not 
in close proximity to the primary excitations. The study was limited to a single point of 
control and should be expanded to multiple actuators and multiple frequencies to evaluate 
the performance of the control. In addition, the control was applied at a frequency close to 
both a structural and acoustic resonance. Other cases of interest include frequencies where 
at most one of the domains is very close to a resonance. 
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Chapter 7 

Conclusions and Next Steps 


A novel active vibration control strategy called Minimum Actuation Power (MAP) is pre- 
sented for a structure excited by multiple primary (excitation) sources and controlled by 
multiple secondary (control) sources. MAP is an optimal control strategy that minimizes 
the total input power into the structure by monitoring the input power from the secondary 
source only and therefore does not need any knowledge about the primary sources. 

The theory of MAP was first developed for a structure excited by a single primary source 
and controlled by a single secondary source. It was shown that the real input (mechanical) 
power in the secondary source goes exactly to zero when the total input power from both 
the primary and secondary source is minimized. This leads to a novel way to implement 
MAP control without explicit knowledge about the magnitude and location of the primary 
excitation. In order to validate the theoretical findings, a model for a simply supported plate 
excited by a primary piezo-electric (PZT) patch actuator and controlled by a secondary piezo- 
electric patch actuator was developed. An electromechanical coupling model of piezoelectric 
patch actuator with a plate structure is developed. A simple relation between the real 
electrical power and the real input power of the piezo-electric actuator was developed. This 
enabled implementing MAP control using the real electric power in the secondary actuator. 
The effectiveness of MAP control was demonstrated for resonant and off-resonant excitation 
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frequencies. It was first shown through simulations that optimizing the location of the 
secondary actuator can help control otherwise uncontrolled modes very effectively. The MAP 
theory was extended to multi-frequency excitation and simulations were presented to validate 
the theory. A test setup that mimics simply supported boundary conditions was built to 
validate the theoretical findings. To implement MAP using electrical power, tests were done 
to characterize the power consumption of the plate and the dielectric loss of the piezo-electric 
patch actuators. The MAP control algorithm was tested for on-resonance and off-resonance 
frequencies and effectiveness of MAP for vibration reduction. We demonstrated up to 47 dB 
vibration reduction at resonance and 7 dB vibration reduction at off-resonance, as measured 
by a single accelerometer at the plate center. There was good qualitative agreement between 
the experimental results and theoretical findings. 

The MAP control approach was extended to multiple primary excitations with zero relative 
phasing and multiple secondary excitations. It was shown that when the primary excitations 
have the same phasing, then the real secondary power is zero when the total input power is 
minimized. However, when the primary excitations have arbitrary phasing, the secondary 
real power is not always zero when the total input power is minimized. It was shown that 
the MAP condition is exactly satisfied when the secondary sources are located symmetrically 
with respect to the primary sources. However, such a constraint on the secondary source 
locations is hard to satisfy in practice. Extending the single actuator model, we developed 
an electro-mechanical model for a simply supported plate excited by multiple primary and 
multiple secondary PZT patch actuators. Performance of MAP control as a function of the 
location and number of the secondary sources was assessed. It was seen that increasing the 
number of secondary sources improves vibration reduction. Also, multiple secondary sources 
can be useful when vibration reduction over a frequency band is sought. Experimental val- 
idation of MAP control was demonstrated in the lab with two primary and two secondary 
excitations using an aluminum plate as well as a composite panel. For example, we demon- 
strated 20 dB vibration reduction for the aluminum plate and 12 dB vibration reduction for 
the composite plate, as measured by a single accelerometer at the plate center. A test for 
real-time implementation was developed to automate MAP control strategy. 

The feasibility of MAP control approach for reducing airframe vibrations due to the rotor- 
craft main rotor transmission and the reduction of corresponding sound levels in the passen- 
ger cabin was assessed using a more meaningful metric for vibration reduction, the global 
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kinetic energy of the plate structure. MAP feasibility was assessed for multiple primary 
excitation sources with arbitrary relative phasing and with multiple frequencies controlled 
by multiple secondary sources. In order to circumvent the problem of placing the secondary 
actuators symmetrically with respect to the primary sources to satisfy the MAP condition, 
it was shown that placing the the secondary sources at the intersection of the node lines 
of the neighboring modes can guarantee that the secondary real power is zero at optimal 
conditions even when the primary sources have arbitrary phasing. An algorithm to compute 
the optimal secondary phase and magnitude was developed for resonant and off-resonant 
excitation frequencies. A stable iterative solution for finding the MAP solution for multiple 
secondary sources was developed. The stable iterative solution, however, is not feasible for 
large number of secondary sources. An alternative feasible approach was developed to find 
the MAP solution when the secondary sources are large in number. A test rig to mimic 
rotorcraft construction was designed, modeled and built. MAP control with primary exci- 
tation using shaker, in addition to PZT, with multiple frequencies and multiple secondary 
actuators for control was demonstrated. It was shown that MAP control achieved significant 
vibration reduction with average global vibration reduction of the order of 15-20 dB for 
single frequency with one secondary actuator. A single secondary actuator was optimized to 
control multiple frequencies, in some cases up to 3 frequencies with 8-10 dB average global 
vibration reduction for each frequency. 

A fully coupled vibro-acoustic finite element model of a larger, more complex and modally 
denser box type structure was used as a quick evaluation of the potential of MAP to reduce 
noise in the passenger cabin due to the main rotor transmission excitation. It was shown 
that MAP control resulted in 6.1 dB reduction in the total kinetic energy of the structure 
while the averaged sound pressure levels (SPL) in the acoustic cavity dropped by 13.8 dB 
at resonant excitation. This preliminary analysis clearly shows the potential of MAP in 
reducing interior noise in a helicopter cabin. 


Next Steps 

Possible future directions include establishing physical system design rules for actuator den- 
sity that leverages the highly distributed actuation opportunities offered by MAP and op- 
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timal actuator positions on a generic panel. It is important to quantify the cost/benefit 
of MAP on an aircraft-level system for both tonal and broadband sources using analysis 
and tests. Finally, there is a need for better understanding of the vibration path and the 
structural contributions to interior noise using a detailed diagnostic survey of vibro-acoustic 
problem and the impact of composite structures and improved gear designs or gear noise 
reduction technologies should be studied to assess the interior noise challenges of the future. 
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Appendix A 


Total input power due to primary and 
secondary patches 


Let us assume that the primary and secondary patch edges are discretized into an equal 
number of elements, say q. Let pi and s t represent the coordinates of the center of the i th 
primary and secondary patch element respectively where i — 1, . . . , q. Let M p and M s be the 
primary and secondary moments acting on each element respectively. Consider the point pi 
on the primary patch. The angular velocity at this point is a sum of the angular velocity 
induced by the the primary and secondary patch moments and can be expressed as 

W K = £ W KK + Y, W K „ (A.l) 

3 3 

where W PiPj is the angular velocity induced at p t due to the moment M p at location pj and 
W PiSj is the angular velocity induced at p t due to the moment M s at location Sj. The angular 
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velocities W PiPj and W PiP . can be expressed as 


PiPj 


W, 


PiPj 


w t 


PiSj 


— PpiPj M p 
= Pp -iSj M s 


(A.2) 


where P PiPj and P PiSj are the cross mobility terms given by 
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A^ S n Pi (Pi)^ (Si' 


PiSj 


E ( 2 I • ( 2 2\\ ^ W/ i r n \*j) ( \ o\ 

(t/ujuj +juj(uj n -uj )) — t~, — — — - (A. 3) 

n=1 K ’ m[(u: 2 n - u> 2 ) 2 + UJ^ 2 } 


where r/ is the structural damping, m is the plate mass, oo is the excitation frequency (assumed 
same for primary and secondary excitation) and oo n is the n th resonant frequency given by 





(A.4) 


where n refers to the mode index ( 77 , 1 , 77 . 2 ), E is the modulus of elasticity, p is the Poison 
ratio and p be the density of the plate material. The function ^ (p) is given by 

#<p) = ^ nos ft) cos ( *&) sin M sin ft) sin ( cos ( *&) 

L'x \ L/x J \ J^y J J^y \ J \ Ly J 

(A.5) 

where (x, 7/) are the coordinates of the point p and is the angle made by the moment axis 
with the x axis. Note that for the primary (or secondary) patch, 8 Pi (or8 Si ) = 0, 7 t/2 , 7 t and 
3 tt /2 for all elements i along edges 1, 2, 3 and 4 respectively. 


The total power input due to the primary patch is given by 


p, = )r<= [m; v Wp \ 


(A. 6) 


Using (A.l) and (A.2) in (A. 6), the total input power due to the primary patch can be 
expressed as 

P p = ^Re ( M* p P pp M p + M*P ps M s ) (A.7) 
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where P pp represents the equivalent angular mobility due to the primary patch given by 


Ppp ~ ^ 


PiPj 


* 3 


and P ps represent the cross angular mobility given by 


Pp s 


* 3 


Similarly, the total power input due to the secondary patch can be shown to be 

P, = ^ Re ( M* s P sp M p + M;P ss Ms) 

where P sp = P ps represents the cross angular mobility and 

Pss = ^ ^ P Si Sj 
i 3 

represents the equivalent angular mobility due to the secondary patch. 


(A.8) 

(A.9) 

(A.10) 

(A.ll) 
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Appendix B 


Total Input Power and Global Energy 


The modal response of the structure can be written as 

w(u,x,y) = 

i 

where (f>i(x,y ) is the mode shape, 

( , 1 fi 

qi {uj) = — — 
m 

is the complex modal amplitude, m is the mass, 

fi = Yl F k(u)(j>i( x k,yk) 


is the modal force and 


^i(w) = w?(l +jr]i) -u 2 


(B.l) 
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is the modal frequency factor where c Oi is the natural frequency and rji is the loss factor. The 
overall time averaged kinetic energy and potential energy of the structure can be written as 
[41] 


E kin = ^J>| 2 (B-2) 

i 

®po. = jE”'. 2 |9‘| 2 (B-3) 

i 

Similarly, the total input power can be expressed in terms of the modal amplitudes as 

— TflUJ V — > 9 9 

p = ~ 2 ~ wi?<i ( B - 4 ) 

i 

Iin general there is no simple relation between P and E^n and E pot . However, if the forcing 
frequency is close to the natural frequency, then the total input power can be relation to the 
total vibrational (kinetic) energy. If u ~ cu;, then a single mode dominates the response. In 
such a case, we have 


Thus, we have 


Thus 


— — mcof . 

Ekin Ep 0 t rs -' ^ 

Qi\ 2 

(B.5) 

P ~ 2 Vi^M 


(B.6) 


P 

(B.7) 

P = 2 TjiUliE fan ^ Eki n = 

2 rjiUJi 

min P = min Eki n 


(B.8) 


That is, minimizing the total input power is equivalent to minimizing the total vibrational 
energy. 

When the excitation frequency is not close to the resonant frequency, that is, co ^ u;*, then 
we can relate the total input power to the potential energy of the structure. Assume that 
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the loss factor is the same for all modes. That is, r)i — rj for all i. In this case, comparing 
(B.3) and (B.4), we get 

P = 2r]ujE pot (B.9) 

The net input power is proportional to the potential energy at all frequencies. Thus, in this 
case 

min P = min E pot (B.10) 


That is, minimizing the net input power is equivalent to minimizing the total potential 
energy 

If there is high modal overlap, it can be argued that 

Ppot ~ Pkin 

for all frequencies except below the first resonance frequency. The argument is as follows 
[41] 

“Below the resonance potential energy of a particular mode is larger than kinetic, while the 
opposite applies above the resonance. In the case of high modal overlap, at any particular 
excitation frequency, the (lower) kinetic energy of modes having the eigenvalues above the 
excitation frequency will become compensated by the (higher) kinetic energy of lower modes. 
Inverse compensation effect will apply to potential energy. As a result, the global kinetic 
and potential energies in modally overlapped systems can be shown to match closely at all 
frequencies except below the first resonance where the compensation cannot take place" 


That is, 


for all u > uq. Thus, 


P m 2r]uE kin 


min P = min E k i n 


(B.ll) 

(B.12) 


Thus, in general, if there is high modal overlap and if the loss factor is same for all modes, 
then minimizing the total input power is equivalent to minimizing the total vibrational 
energy for all frequencies except below the first resonance frequency. 
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